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1. Introduction 


Applied Research, Inc. is pleased to submit the enclosed 
final report which concludes contractual efforts on Research 
Study: Beta Experiment under NASA/MSFC Contract NAS8-3^337. 

The contractual effort had a period of performance extending 
to July 15, 1982. 

This contract included twelve (12) tasks related to the 
design, development, calibration and test of a focused Laser 
Doppler Velocimeter (LDV) system to be employed for the 
measurement of atmospheric backscatter (Beta) from aerosols 
at infrared wavelengths. 

Task 1 involved the design, development, and fabrication of 
a Doppler signal generator which were used in mapping the 
coherent sensitive focal volume of a focused LDV system. 

Task 2 provided technical assistance for design and 
fabrication of optical components that were employed in the 
Beta measurement system. 

Task 3 provided technical assistance for the analysis of 
system calibration data during the flight test activity 
scheduled for the Beta system during July 1981. These 
analyses were performed to determine the acceptability of 
the Beta measurement system's performance and to make 
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recommendations on the operating parameter setting to 
increase system performance when needed. 

Task ^ provided for support during the system integration 
phase into the NASA/ARC CV-990 aircraft. This support 
included ground testing and system calibration after the 
equipment was loaded into the aircraft but before the start 
of the flight test series. 

Task 5 involved the design of a Beta measurement system 
which operates at 10.6 /xm and 9.1 Mm simultaneously. The 
technical effort included design trade-offs, hardware and 
software recommendations, and an assessment of the degree to 
which common hardware could be shared by the two wavelength 
systems . 

Three contract modifications were attached to the NAS8-3^337 
contract. These modifications included seven additional 
technical tasks to be performed within the existing period 
of performance. These modifications included seven 

additional technical tasks to be performed within the 

( 

existing period of performance. These tasks were referred 
to as Tasks 6, Task 7, 8, 9, 10, 11 and Tasks 12. 

Task 6 involved the providing of manpower, facilities, and 
materials to accomplish the evaluation of the applicability 
of adaptive filtering to signal processing of the data from 
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NASA/MSFC pulsed Doppler lidar. This task specifically 
analyzed the applicability of adaptive filtering to Doppler 
lidar data from calibrated targets, aerosol targets with 
known velocity, and aerosol targets with turbulence. 

Task 7 involved the development of software algorithms, 
experiments, data processing capability, and data analysis 
capability as required for the correlation of the Beta 
Measurement data with previously recorded laser Doppler 
system data and related aerosol density data. 

Task 8 involved the design and performance of a sequence of 
measurements to determine Lgff (the effective range 
increment to be associated with a hard target signal), the 
determination of the ranges at which these measurements 
should be made, and the determination of the supporting data 
(such as frequency of Doppler shift and lidar system output 
power). These data were analyzed and the results documented 
in regular monthly reports. 

Task 9 involved the design and performance of a test which 
utilized a wire calibration target which was spun through 
the lidar beam. The detected signal was digitized and 
processed to provide sensitivity contours which were 
employed in the calibration of the Beta system. The results 
were reported in regular monthly progress reports. 
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Task 10 involved the defining of a technique by which the 
response of the lidar to a single aerosol could be measured. 
This included design and fabrication of apparatus for 
generating and propelling single aerosols through the lidar 
beam. Development of software to analyze the resulting 
lidar data, along with documentation of the measurement 
results were included in this task. 

Task 11 involved development of software which was required 
to transfer backscatter data collected during the summer 
1981 flight program to the Sigma V computer. The 
backscatter data which existed in protocol, file structure 
device characteristics, signal levels, etc., at NASA/MSFC 
were provided GFE to Applied Research, Inc., personnel. 

Task 12 involved the generation of the necessary software to 
interface a Hewlett-Packardd lOOOL computer with several 
peripheral devices for the purpose of transferring, 
processing, displaying, and storing of recorded data. 
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2. Task 1: Doppler Signal Generator for Coherent Focal 


Volume Mapping 

ARI has designed and assisted NASA/MSFC personnel in the 
development and fabrication of a Doppler signal generator 
for use in mapping the coherent sensitive focal volume of a 
focused LDV system. This apparatus consists of a method of 
sending single aerosol particles into the focal volume 
region, called a single particle generator, and a spinning 
wire target. Figure 1.1 shows a photograph of the completed 
single particle generator. 

The single particle generator consists of a plexiglass test 
chamber, a pneumatic control panel, a particle generating 
chamber and a particle injector. These units were 
constructed early in 1981, assembled and tested in 
September. These preliminary tests revealed problems in the 
design that required modification to correct. 

1. The Test Chamber: The test chamber is constructed of 

5/16" acrylic sheet. This was originally assembled with 

epoxy. However, under negative pressure the long seams 

proved to be too weak, resulting in leaks. This was 

corrected by reassembling the structure with 1/2" aluminum 

reinforcing members at all seams. The plexiglass is 

fastened to the aluminum bars with 5-32 screws and all seams 

are sealed with RTV. A new top was also constructed from 
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1/2" acrylic sheet. Since the top must be removable the 
added thickness was required for rigidity. Aluminum bars 
were fastened to the top to support the inside top edges of 
the sides. The wing nut fasteners for the top were 
originally secured to the sides of the chamber by epoxy. 

When negative pressure was applied to the chamber the sides 
flex enough to break the epoxy bond. This was corrected by 
using screws in addition to the epoxy to secure the 
fasteners. The flexing of the sides also caused the gasket 
providing the seal for the germainum window to leak. A 
channel was milled around the window opening and the gasket 
was replaced by an o-ring. 

2. Pneumatic Control Panel: The pneumatic control panel 

was originally configured as shown in Figure 1.2. This 
configuration was based on the requirement for the system to 
operate at either positive or negative pressure. Since the 
system will now only operate at negative pressure and 
because of modifications to other parts of the system, the 
control panel has been reconfigured as shown in Figure 1.3. 

A Millipore filter has been added to the air inlet to reduce 
the number of ambient particles in the test chamber. 

3. Particle Generating Chamber and Injector: Originally 

the particle injector was a commercial air brush connected 
by tygon tubing to the generating chamber. Although this 
had worked well in tests before the test chamber was 
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constructed, it did not function properly when placed in the 
negative pressure environment. The air brush was replaced 
by a hypodermic needle attached directly to the tubing. 

This arrangement worked, but not reliably. There was also a 
considerable build up of particles in the tubing which 
occasionally clogged the line. The generating chamber was 
moved inside the test chamber and the hypodermic needle was 
fitted directly through its side, giving direct access to 
the particles (see Figure 1-4). The generating chamber is 
an aluminum cylinder with a speaker in the bottom to provide 
sonic excitation (at about 600 Hz) to the particles. Aside 
from the relocation and modification mentioned above it is 
unchanged . 

In addition to these modifications, a water manometer was 
built and connected in the test chamber vacuum line. The 
manometer permits more precise measurement of the pressure 
in the test chamber and thus more control over the air flow 
through the particle generating chamber. The exit velocity 
of particles from the single particle generator may be 
estimated from Bernoulli's equation. The following figure 
shows the assumed model in which the density of the fluid is 
taken to be constant. 
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Because of the difference in areas at points 1 and 2, Vj is con- 
sidered negligible compared to V 2 . Hence 
5 v/ - 2 P / 

P . = .081 lb / ft^ 

air 



= 1.51 X 10^ h 

V = 19.6 Vp (m/sec, with h in inches) 

Under measurement conditions h was found to be about .125 
inch. Hence V = 6.9 m/sec is the expected velocity, under 
the above assumptions. 
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3. Task 2: Assistance for Design and F abr ication of Lens 

Mounts 

Technical assistance in the design and fabrication of 
optical components was provided in the fabrication of a 
gimbal mirror mount for the output mirror of the telescope. 
Modification of the case which encloses the interferometer 
and laser was accomplished in order to provide improved 
access to controls. Optical components including lens 
mounts, detector mount fixture, zinc selenide (ZnSe) lens 
and adjustable lens mount fixtures were designed, fabricated 
and procured. The Beta system design has the support 
structures to accommodate a flexible design which uses 
either the Honeywell or the Hughes CO 2 laser. 

A decision to employ the Honeywell CO 2 laser for the initial 
flight test series permitted optical support fixtures for 
optical elements to be fully integrated into the Beta 
system. The Beta system was assembled and preparation was 
made for calibration tests, prior to shipment, for the 
volume backscatter mode. 
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4. Task 3: System Calibration Data Analysis 


During July, 1981, Applied Research, Inc. personnel provided 
technical assistance for the analysis of system calibration 
data during flight test activities onboard the NASA/ARC 
Convair 990 aircraft. These analyses were performed to 
determine the acceptability of the Beta measurement system’s 
performance and to make recommendations on the operating 
parameter settings to increase system performance. 

To determine /8 in the volume mode, the expression S = )3 G is 
used. For an untruncated Gaussian beam it has been shown 
(T.R. Lawrence et .£l , RSI , 512 ( 1 972 ) ) that G =-rv\?/B, 

where 

v= system efficiency 
A = wavelength 

P = number of transmitted photons/sec 
B = electronic bandwidth. . 

For truncated beams an appropriate system response must be 
used. This response can be represented by the same 
expression, except that v becomes a function of focal range. 
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The volume mode calibration involves determining the system 
response G to a volume distribution of aerosols for a beam 
which is usually significantly truncated. In principle, G 
could be determined by mapping the focal volume with a 
single particle probe and integrating. In practice, 
calibration is usually accomplished by observing the 
backscattered signal from a plane rotating disk. For the 
truncated Gaussian beam, the volume signal S and disk signal 
Sj) are related by 

Sj) = AL S with aL = f^ X /R^ 

where f is the focal length and R is the Gaussian radius. 
From this it follows that p= AljS, where p is the 
bidirectional reflectance of the disk in inverse steradians. 
The length A L is a range increment from which approximately 
50% of the energy is returned in the volume mode for an 
untruncated Gaussian beam. Hence for such a beam, the 
efficiency factor may be calculated from the length AL and 
the known reflectance of the disk. To calibrate a truncated 
beam the corresponding aL must be found by measurement or by 
computation of truncated beam propagation characteristics. 
This quantity is found to vary differently with range from 
the untruncated AL. 

Volume mode calibrations of the Beta system were 
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accomplished using both the sandpaper and sulfur Doppler 
targets. The system performance from the laboratory test 
showed 68 dB si gnal-to-noise from the sandpaper target and 
82 dE signal-to-noise from the sulfur target. These signal- 
to-noise figures indicate an overall system efficiency of 
approximately 105 for the Beta system operating with 6 watts 
transmitter power and a spectrum analyzer bandwidth of 100 
KHz. These performance numbers indicate that the current 
Beta system will have sensitivity to measure atmospheric 
backscatter to approximately 10"^® m"”* str“^ with a signal- 
to-noise of 1. 

In June, 1981, the Beta system was used to observe the 
atmosphere at a focal range of approximately 100 ft and with 
line of sight winds greater than 2 m/sec. SNR's in excess 
of 20dB were observed with a transmitted output power of 6.2 
watts . 
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5. 


Task Assistance in Flight Preparations and Ground 
Testing 


Applied Research, Inc. personnel assembled the protective 
cover for the Beta system which will be employed to meet air 
safety requirements when flying on the CV-990. The Beta 
system with all supporting hardware and equipment was 
packaged and shipped to ARC for participation in the Severe 
Storms Measurement System flight test program. 

Work continued at Ames Research Center (ARC) to provide 
technical support during the system's integration phase into 
the NASA/ARC CV - 990 aircraft. This support included 
ground testing and system calibration following shipment of 
the system to ARC before CV - 990 flight test activities 
began . 

ARI assisted in fitting of the equipment into an aircraft 
flight rack. This included the spectrum analyzer oscilloscope, 
chiller, stepping motor monitor panel, CRT terminal, signal 
processor, disc drive and laser power supplies. The stepping 
motor junction box was interconnected with new cable. ARI also 
assisted in testing the system by measuring and analyzing the 
return utilizing the shutter meachnism. 
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6. Task 5: DESIGN OF A TWO-COLOR BETA MEASUREMENT SYSTEM 


1. IKTROD'JCTION 

The recent success of NASA/KSFC's Laser Doppler Velocimeter System, 
which is designed to measure atmospheric backscatter, £, at 10.6 ym has 
generated an interest in determining the design of a two-color 6 measure- 
ment system. Applied Research, Inc. has undertaken an analytical system 
design effort, centered around existing 6-system hardware, to result in 
a two-color 6 system. This effort is focused toward employing e>i sting 
hardware to the maximum extent possible and does not explore the develop- 
ment of new signal processing electronics and/or mixing of two-color 
reflected signals on a common detector with two local oscillator frequencies. 
The results of this design study, including a ray trace analysis of a 
proposed two-color 6-system design which employs the current system's beam 
expander, is given in the following sections of this report. 

2. DESIGN’ CONSIDERATIONS 

The decision to investigate two-color -system designs that employ 
the maximum amount of existing MSFC hardware has been driven by consideration 
of other efforts which are more extensive and will explore multi-color 
-system designs. The design to be derived from this study was not 
considered to fit inside a pre-detemined small volume envelope or weight limit. 
This design effort considered several two-color system designs which \ilti- 
matelj’ result in two design options: 
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• A bea^ splitter combined two-color system which 
results in a common focal volume for the two transmitted 
beams. (System A) 

• A two beam directing insertion mirror system which 

uses a common beam expander and focuses the two transmitted 
beams at approximately the same spatial location. (System B) 

System A was analyzed for two configurations. One design which em- 
ployed a 50/.' beam splitter as the beam combining component (for and 
X 2 ) on the transmit leg of the system was analyzed and found to have 
relatively low performance potential when compared to the other designs. 

The other design which was considered employed a dichroic beam 
splitter as the beam combining component. This design was analyzed and 
found to have relatively higher potential for performance than the 50% 
beam splitter design. Other considerations employed a dichroic beam 
splitter in the return leg to separate reflected signals (X| and X^) 
such that they could be combined with separate LOs and heterodyned on 
separate detectors. These designs are more fully discussed in Section 3.1, 
Design Analysis of the Beam Splitter Combined Two-Color 6-Systems. 

The second design, Two-Beam Directing Insertion Mirror System, 

System B was analyzed and found to have a very good system performance 
potential relative to the existing single-color 6-systems. This second 
design was chosen for further system design analysis because it can poten- 
tially employ the existing 6-system‘s interferometer, laser transmitter, 
beam expander, detector, and signal processing. This design will require 
the use of an additional signal processor (or time-sharing with the current 
signal processor), laser of wavelength X 2 » detector, interferometer and a 



folding/bear, directing flat mirror for beam insertion of the second 
color onto the existing beam expanding secondary. This system's design 
will be such that the spatial relationship of the two focal volumes 
will be known for various focal ranges and such that statistical correlation 
of s(r) may be made for the two-colors whether employed in the single 
particle mode or the volume mode. 

The System B design obviously employs more hardware and thus has larger 
volume and weight characteristics than the System A design. However, from 
the design analysis which predicts system performance, it is indicated 
that the System B design will provide greater sensitivity (fewer system 
losses) assuming that all other system characteristics such as laser power, 
optical efficiency, heterodyne efficiency are held constant. 

3. DESIGN ANALYSIS 

Several two-color B-system designs were analyzed to determine their 
performance characteristics. All the designs were found to reduce to two 
basic two-color B-system designs which will be analyzed in Sections 3.1 
and 3.2. 

A 6-system which employs a beam splitter for combining the two 
transmitter output signals is analyzed in Section 3.1. A B-system which 
employs two folding insertion mirrors is analyzed in Section 3.2. As a 
system design trade-off between the two potential designs, a ray trace 
analysis which incorporates the effects of off-axis insertion of one of the 
input beams, is performed. The results of this analysis indicate that 
the signal to noise loss due to off-axis insertion of the second color is 
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less expensive, in terms of s/n loss, than the 50 % beam splitter com- 
bining system: which has 3 db of losses on the transmit leg and approximately 
3 db of losses on the receiving leg of the system. 

The dichroic beam splitter version of the System A has potential for 
an improved s/n performance over the 50« beam splitter version. Initial 
discussions^^\vith vendors of such beam splitters indicate that as great 
as 80/i transmission of 10.6 radiation can be achieved while at the 
same time reflecting as much as 80% of the 9.1 ^m radiation. This per- 
formance v.'ill improve the System A design performance significantly both in 
the transmit leg and in the receiver leg, through use of separate L.O.'s. 

The transmitter and receiver leg loss are each reduced to less than 1 db. 

3.1 DESIGN' ANALYSIS OF THE BEAM SPLITTER COMBINED TV/O-COLOR 

Two options are presented with respect to beam splitter choices for 
this design. The first option employs conventional beam splitters, while 
the second option utilizes dichroic beam splitters, at critical points, 
in order to enhance the transmittance and reflectance at Xj and 
desired. 

3.1.1 CONVENTIONAL BEAM SPLITTERS 

Fig. 5.3.1a gives the optical configuration for the beam splitter 
combined two-color 6-system. This system combines the radiation from 
two CO2 Lasers (X^ and Xg) and passes the combined beam through a 
single interferometer. The output beam from Laser 1 is reflected off 
mirror and onto beam splitter Sj. The radiation from Laser 2 is re- 
flected on beam splitter such that the output powers of the two lasers 
are combined at Sj and are of approximately equal power. The requirement 
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to have equal power components of >.j and >.2 in the combined transmitted 
beam results in approximately 3 db of output power loss from each of 
the laser transmitters. The radiation that is reflected by must be 
dumped from Laser 1 and the radiation from Laser 2 that is transmitted 
by Sj must be dumped from the system. Fig. 5. 3. lo shows this radiation 
being dumped into a radiation sink. The remainder of the system is exact- 
ly as the 6-system which flew on the Ames Research Center's Convair 990 
as part of MSFC Summer 1981 test program. (Reference Fig. 2, "Design and 
Calibration of a Coherent Lidar for Measurement of Atmospheric Backscatter", 
SPIE Conference, San Diego, CA. May 1981.) A positive attribute of this 
system is that the two colors are focused in the same spatial volume. This 
permits comparison of single particle signal returns for single particle 
measurement operations and a direct comparison of backscatter for the volume 
mode operation (without applying statistics to the signals). 

This system may be operated in one of two possible modes for handling 
the signal processing of the reflected returns. The reflected radiation 
Xj and X 2 may be heterodyned with the combined beams local oscillator signal 
which will result in the use of one detector and signal processor. The 
Doppler returns when mixed with the combined Xj and Xg L. 0. radiation 
can be sufficiently separated such that the signal processing of the separate 
colors will be possible. However, the L.O. radiation of X^ will appear 
as noise for the processing of X 2 (reflected returns for X£) radiation 
against the L.O. and vice-versa for and X^* This effect will tend to 
reduce the operating efficiency of this design further, maybe as much as 3 db. 
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This r-system design suffers from a system efficiency point of 
view in that the inherent losses for each color may be as much as 6 db, 

3 db on the transmit leg and 3 db on the receiver leg. 

3.1.2 DICHROIC BEAM SPLIHERS 

A dichroic beam splitter can be designei^^ which will permit high 
reflectance at transmittance at Xj, or vice-versa. In Fig. 5.3.1b 

such a splitter is utilized to effect separate mixing of Xj and 
separate detectors. 

Transmit beams from Laser 1 and Laser 2 are combined at $ 2 . a dichoric 
splitter with strong transmittance at and strong reflectance at X 2 * 

Local oscillator signals from Laser 1 and Laser 2 are not mixed, as in 
Fig. 5. 3. la, but are combined separately with returned signals at and 
Sg, which are dichroic splitters with high reflectance for and ^2 
respectively. It is only necessary to obtain enough L.O. upon transmittance 
to maintain shot noise limited operation. Also, the effect of strong re- 
flectance of the Xj L.O. onto Sg will be minimized since Sg has strong Xj 
transmittance. This system offers about 2 db less power loss in transmit 
and receive legs than the previously discussed option. 

3.2 DESIGN ANALYSIS OF THE 6-SY$TEH EMPLOYING TWO BEAMS WITH FOCUSING 

AT SEPARATE POINTS 

The schematic design of the two-beam system is shown in Fig. 5.3.2. 

This design has inherently good optical efficiency in that it avoids the 
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FIGURE 5.3.2 

TWO COLOR BETA SYSTEM OFF-AXIS BEAM INSERTION DESIGN 


previously discussed system's 3 db output power loss on the transmitter 
ooticel leg and avoids the 3 db signal loss in the receiver optical leg 
of the system. However, the potential for heterodyne efficiency degradation 
due to off-axis insertion of the X 2 beam must be investigated to insure 
that the design's overall system efficiency is preserved. The following 
is a design analysis of the two beam directing insertion mirror system 
which considers the proposed optical design and employs the existing 6-system 
optical beam expanding telescope. Operating ranges of 10 meters and 30 meters 
have been considered in this analysis. It should be noted that focussing of 
this system at longer ranges tends to reduce the effect of off-axis beam in- 
sertion. 

As depicted in Fig. 5. 3. 2 this optical design results in the two colors 
being focused at approximately the same range but in slightly separated 
spatial volumes. 

A detailed discussion of this system's optical design and ray trace 
analysis is as follows. 

One of the two flat mirrors, Mj and shown in Fig. 5. 3. 2, is 
located above the optic axis (out of the paper), with Xj parallel to 
the axis. The second X 2 beam is located in the same plane (parallel 
to the paper) but displaced as shown. The plane of the two laser beams 
striking the secondary (small) mirror, is perpendicular to the plane 
of bilateral symmetry of the beam expander. The angle between the two 
beams is made as small as practicable to minimize aberation of the X 2 
beam. 
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The center of the beam at the flat mirror is 1.3 cm from 
the plane of bilateral symmetry through the optic axis (perpendi- 
cular to the paper). The mirror is then adjusted so that both 
beams have spatially identical images on the large focusing mirror, 
(called the primary). These images are approximately 4 inches 
in diameter, as shown in Fig. 5. 3. 2. Note that the flat mirrors 
appear to obscure the beams coming from the secondary (see figure) 
but in fact do not, since the primary is completely off axis, in a 
direction out of the paper. 

The center of the beam strikes the secondary slightly above 
the optic axis (out of the paper). The center of the beam strikes 
the secondary at the same distance above the optic axis and at a 
point that is displaced slightly in a perpendicular direction 
(toward the top of the paper). The amount of this lateral displac- 
ment of the point on the secondary is, however, less than 0.1 mm, 
and can easily be ignored in ray-tracing calculations. To see this, 
note that if the central ray of X 2 struck the secondary at the same 
point as that of X^, it would be displaced at the primary by about 
1.4 mm from the central ray of Xj. 

The Gaussian focal point (aberrations ignored) of the off- 
axis beam will be in the focal plane of the on-axis beam, and can 
be obtained by a simple geometrical construction. It is displaced 
to the opposite side of th6 optic axis from the plane mirror, and 
the displacement distance for the 10 meter focus is 0.3 inches. 

By comparison of the wavefront shift at different points on 
the secondary in Tab. 5. 3. S and Tab. 5.3.10, we see that the two beams 
should have essentially the same character! sites. A second beam 
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added to the Lidar system, as described in this section, should give 
results almost as good as the original beam. Indeed, the differences 
between the two would be inconsequential. 


3.2.1 ABERRATIONS OF THE TVJO-BEAM SYSTEM 

The geometrical theory of aberrations is given by many authors. 

The authoritative source Born and Wolf ^ (B & W) was used for the 
calculations reported here. The important performance characteristics 
of an optical system may be more adequately analyzed for the purpose of 
this design effort using the Born & Wolf approach. Most of the pertinent 
material is contained in Chapter V, pp 203-232. References to pages numbered 

100 + are to pages in B and W. There are five third-order aberrations 
of a monochromatic system: Spherical aberration, astigmatism, curvature 

of field, distortion, and coma. Associated with these are five co- 
efficients - B, C, D, E, and F, respectively - which can be calculated 
from the Gaussian variables (object distances, radii of curvature, image 
distances, etc.). 

For a two mirror system like the beam expander considered here, 
the object distances are defined to be and Sg; image distances, Sj and 
S^; and radii of curvature, r^ and With the proper choice of signs 
for radii of curvature, (both positive) and indices of refraction 
(n^ = 1, nj = -1, Og ® 1), one gets the following related quantities 
(p. 224): 

h, = -1 

K = -SJS\ 

kf = 0 ^ ^ 
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where d, is the di stance between the primary and secondary mirrors. 
Four other quantities involved in the calculation of the aberration 
coefficients are 


K = 1 .1 

1 r, Sj 

K . i .i. 

2 =2 


Eq. (15), p. 223 


and b^ = -1 = b^ The values of bj and b^, the deformation coefficients 
for the mirrors, represent the fact that the mirrors are paraboloidal. 

It is convenient to define three more sets of intermediate variables 

“i = 1 *>i ("i - 

Y, = ! ( l/n,sV - Si')- 

The equations for the aberration coefficients (Eq. (24), p. 225) then 
become 

4 2 

B = I [ a.+ h . K.T. ] 

C = Z [ a. + ( 1 + B. )^ T. ] 

D = T [a. + 6^ (2 + B. ) Y. ] 

E = ^ [a . k? +k. ( 1 + B. ) (2 +B. )Y. ] 

F = T fo,. k^. 4h2 K. ( 1 +B. )Y . ] 
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The aberrations themselves are calculated using the so-called Seidel 
variables (pp. 20S), For the purposes of .this design analysis these are 

^0 = ''c / 

where is the lateral displacement of the object point from the optic 
axis and is the object distance to the first mirror, and 
^ 2 = {^2 + D 2 P 2 )/M' 
n 2 = (Y 2 + D 2 d2 )/M‘ 

where (X 2 , Y 2 ) is the lateral position of the image point, (-D 2 ) 
is the image distance from the second mirror, and M' is the lateral 
magnification between the two mirror planes. The ray components P 2 and q 2 
are direction cosines at the image plane (p. 134) with respect to the 
X and y axes, respectively, where the optic axis is the 2 axis. 

The displacement of the wavefront (fourth-order correction) and 
displacement of the image point (third-order correction) are then cal- 
culated. Let ^2 ” 3nd Yj - Yj* be the corrections to the position 

of the image point in the image plane ( p.205). 

Define p and 6 by 


^ 2 ~ P sin0 (p. 213) 
n 2 * P COS0. 

and define the displacement of the wavefront due to the aberrations 
as . Then with 




(3) 


X 



one finds (Eq. (7-8), p. 213) 
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< - i - C cos - i Dy^ + E y^ pcosl 

3 

+Fy^ p cos3 

3 2 2 

' X = B p sin G-2 Fy^ p sin G cos G + D y^ p sin G 

y = B p^ cos G- Fy^ p^ (1+2 cos^ G) + (2C + D) y^ pcosG - E 

thus, there exists a complete method of calculating the aberrations 
once the Gaussian parameters are known. 

3.2.2 BASIC PARAMETERS OF THE TWO-BEAK SYSTEM 

The basic parameters needed for geometrical calculations are for the 

radii of curvature of the two mirrors and the distance between them. The 
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basic data were given directly from the blueprint ' supplied by the 
original manufacturer of the beam expander. The most accurate data are 
the original distance between mirrors (for parallel beams) and the magni- 
fication which are respectively 32.670 and 14.987 inches. [Note: Hence- 
forth, unless otherwise noted, all lengths are in inches. To avoid any 
round-off errors, more than the actual number of significant figures 
are given for some numbers.] The focal lengths are then determined to be 

fp (primary) = 35.006 
fg (secondary) = 2.3357 
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Most of the remaining parameters depend on the location of the image 
plane. The design effort considered examples where the image planes are 10 
meters and 30 meters from the primary mirror. These are actually the focal 
planes of the two-mirror system. 

3.2.3 ABERRATIONS IN THE FOCUS OF THE INITIAL LASER BEAMS 

There are two basic problems to consider: (1) The aberrations 
in focussing the original beams and (2) the aberrations in the return- 
ing beams. The first case is considered in this section. The distances 
of the image plane (focal plane) from the primary (large) mirror are 10 
meters (393.7 inches) and 30 meters (1811.1 inches). 

The quantities needed for calculating the aberration coefficients 
(Section 3.2.1) are given in Tab. 5.3.1. 
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TABLE 5.3.1 


Quantities Needed for Calculating Aberration Coefficients 
for 10 - meter and 30 - meter Focus . 


a) 

Quantity 

10 Meters 

30 Meters 

b) 


’'l 

4.6715 

4.6715 

'~2 

70.012 

70.012 

^1 

36.087 

33.740 


-00 

-00 

/ 

2.3357 

2.3357 

^2 

38.423 

36.076 

^2 

393.7 

1181.1 

M' 

16.46 

15.44 


a) See Section 3.2.1. 

b) All lengths in inches. 


From these quantities one can determine the related Quantities hp 
kp Kp Qp 6^ and the aberration coefficients (see Section 3.2.1). For 
the coefficients one gets the results in Tab. 5.3.2 
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TABLE 5.3.2 


Aberration Coefficients for Focus 
of Initial Laser Beams 


a) 

Coefficient 10 Meters . 30 Meters 

B -0.069152 -0.019091 

C -0.73171 -0.49090 

D -0.53193 -0.29112 

E 1.3227 0.79201 

F 0.15127 0.041664 


a) See Section 3.2.1. The dimensions of the different coefficients 
are not uniform. Numbers given are in units of (inches)*^ where 
n is an integer. 


Rays to the image point can be identified in terms of points 
on the primary mirror from which they were reflected. This is the 
easiest way to determine the ray exponents Pg and ^ 2 * which are 
the direction cosines (Gaussian approximation) of the rays at the 
image plane. To encompass the beam, we pick points around the 4 inch 
spot made by the beam on the primary. Let the y - z plane be the plane 
of bilateral symmetry of the beam expander. We then label rays in terms 
of X and y positions on the primary mirror as follows: 



X Position 


Y Position 


Ray 1 

5.80 inches 

0 

Ray 2 

3.80 inches 

2.00 inches 

Ray 3 

1.80 inches 

0 


These labels will be used in the calculations to follow. 

The Seidel variables are calculated from the equations in Section 3.2.1. 
The angles of the rays are small enough that sines may be replaced by ‘ 
tangents to 3 significant figures. The image points for the Xj beam 
are, of course, on the optic axis. For the beam, it is sufficient 
to find one ray through the image plane to locate the image point in the 
Gaussian approximation. The lateral magnification M' between the two 
mirror planes depends on d^. Its values are given in Tab. 5.3.1. This 
means, of course, that the extent of the beam striking the primary varies 
slightly with image distance, by less than 1 %. This variation has not 

been taken into account in choosing Rays 1-3, but fixed points on the 
primary have been used to define tne rays, for simplicity. 

For the A^ beam, we see that y^ = 0, so only B is involved in the 
aberrations (spherical aberration). For the beam, y^ is the tangent 
of the angle made by a ray in the initial beam with the plane of bilateral 
symmetry of the beam expander. This involves the location of the flat 
mirror, which is taken as 3.00 inches forward of the primary, with the 
beam center striking the flat mirror at 1.3 cm from the symmetry plane. 

This gives 

y^ = 1.31 X 10”^. 

and for the positions of the Gaussian focal points of the oeatn, 
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(X^, Yj) = (0, -0.314) for 10 meters. 

= (0, -1.00) for 30 meters. 

Using these values, the Seidel variables € 2 * ^ 2 * P. sin £, and cos G, 

(4) 

and from them the wavefront shifts (J’ ' and displacements 

* * 

(Xj - Xj, Yj - Yj) can be calculated. The results for the Xj 
beam is given in Tab. 5. 3. 3-5. 3. 4; those for the ^2 
Tab. 5. 3. 5-5. 3. 6. 
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TABLE 5.3.3 


Seidel Variables and Aberration 
Results for On-Axis Beam; 10 - Meter Focus. 


Quantity 

Ray 1 

Ray 2 

Ray 3 

^2 

0.352 

0.231 

0.109 

r,2 

0 

0.122 

0 

c(4) 

6.7 ym 

2.0 ym 

0.1 ym 

X2'X'2 

- 0.072 

-0.026 

-0.002 

★ 

Y - Y 
^2 2 

0 

-0.014 

0 


a) All lengths in inches except as noted. 
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TABLE 5.3.4 

Seidel Variables and Aberration 
Results for On-Axis Beam, 30 - Meter Focus. 


Quantity 


4 



Ray 1 

Ray 2 

Ray 3 

0.376 

0.246 

0.117 

0 

0.130 

0 

2.4 urn 

0.7 urn 

0.03 • 

-0.078 

-0.028 

-0.002 

0 

-0.015 

0 
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TABLE 5.3.5 

Seidel Variables and Aberration 
Results for Off-Axis Beam; 10 - Meter Focus 


Quantity 

Ray 1 

Ray 2 

Ray 3 


0.352 

0.231 

0.109 

r>2 

0 

0.122 

0 

V 

6.7 urn 

2 . 6ym 

0.008 

- X* 
★ 

- 0.072 

-0.029 

-0.002 

Y - Y 

'2 '2 

- 0.006 

-0.018 

-0.001 
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TABLE 5.3.6 


Seidel Variables and Aberration 
Results for Off-Axis Beam; 30 - Meter Focus 


Quantity 

Ray 1 

Ray 2 

^2 

0.376 

0.246 

r; p 

0 

0.130 

(4) 

C 

2.5 ym 

1.0 ym 

* 

Xp - Xp 

- 0.079 

-0.031 

★ 


- 0.006 

-0.022 


Ray 3 
0.117 
0 

0.03 ym 
-0.003 


- 0.001 



The numbers in Tab. 5. 3. 3-5. 3. 6 can be used to define the spot 
size of the focus in each case, as well as the wavefront variation 


for rays coming to different points on this spot. By comparing Tab. 5.3.3 
with Tab. 5.3.5 and Tab. 5.3.4 withTab. 5.3.6, it can be seen that the 
difference brought about by having the Xg beam offset from the synrietry 
plane is indeed slight. In regard to focussing properties, the two- 
beam 6-system should perform essentially as well as the one-beam 
system. 


3.2.4 ABERRATIONS IN THE RETURNING BEAMS 

The question of the effect on waves reflected from the aerosol 
particles in the focus of a beam will now be considered. Again the 
cases of >.j and beams at 10 and 30 meters are treated. The roles 
of primary and secondary mirrors are now reversed. Quantities correspond- 
ing to those in Tab. 5. 3.1 are given in Table 5.3.7. 
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Quantity 

"l 




table 5.3.7 

Quantities Needed for Calculating 
Aberration Coefficients for Returning Beams 
10 Meters 
70.012 
4.6715 
36.087 
393.7 
38.423 
2.3357 
• 00 

0.06075 


30 Meters 
70.012 
4.6715 
33.740 
1181.1 
36.076 
2.3357 

oc 

0.06477 
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The related quantities are calculated as before, and we arrive 


at the aberration coefficients given in Tab. 5.3.8 


TABLE 5.3.8 


Aberration Coefficients for Returning Beams 


Coefficient 

B 

C 

D 

E 

F 


10 Meters 
9.4423X10'^ 
-1.269X10’^ 
-2.0105X10’^ 
-7.5126X10^ 
1.61X10"^ 


30 Meters 
3.3522X10’^ 
-4.21X10’^ 
-2.0020X10"^ 
-6.6612X10"^ 
1.7X10"^ 


f4l 

The parameter, 4' , the wavefront displacement at different 

points on the secondary mirror is calculated from 


^2 = X2' / M‘ 
ri2 = Y 2 / M ' 


where (X^ , is the point at which a ray leaves the secondary. For 
simplicity Ray n' is defined as the reverse of ray n. Results are given 



TABLE 5.3.9 


Quantity 


10 meters 



30 meters 



Seidel Variables and Aberration 
Results for Returning Beam on Axis 


Ray l' 

Ray 2 ' 

Ray 3’ 

5.80 

3.80 

1.80 

0 

2.00 

0 

-6.8 ym 

-2.0 ym 

-0.06 ym 

5.80 

3.80 

1.80 

0 

2.00 

0 

-2.4 ym 

-0.7 ym 

-0.02 ym 
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TABLE 5.3.10 


Seidel Variables and Aberration 
Results for Returning Beams Off-Axis 


Quantity 

Ray l‘ 

Ray 2 ' 

Ray 3 


h 

5.80 

3.80 

1.80 

~10 meters < 

r,2 

-1.51 

2.47 

-1.51 



-7.7 vim 

-2.5 ym 

-0.2 ym 


h 

5.80 

3.80 

1.80 

— 

'•2 

c<^) 

-0.442 

2.44 

-0.442 

30 meters | 

-2.4 ym 

-0.08 ym 

-0.02 y m 
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3.3 HETERODY!;- EFFICIEMCIES FOR THE ON-AXIS AND OFF-A);i$ REFLECTED 


RADIATION : The heterodyne efficiencies associated 

with radiation backscattered into the system from Xj and 
focal volumes can be calculated using the information contained 
in Table IX and X and Fig. 5. 3. 3. Fig. 5.3.3 is taken directly 
from "Heterodyne detection: phase front alignment, beam spot 
size, and detector uniformity", Steven C. Cohen, August, 1975/ 

Vol. 14, No. 8/ Applied Optics, page 1957. The assumptions 
for use of the Figure 5.3. 3 in calculating the heterodyning 
efficiency are as follows: 

0 The local oscillator is assumed to have a Gaussian 
beam profile. 

• The returns signal beam profile is assumed to be an 
Airy Disk. 

c The Airy Received signal parameter is given by 
irr 

yp where r^ (beam radius on the detector) _ .05 mm, a equal 
the wavelength of the transmitted radiation, and F equal 
the F/No. for the focusing lens placed immediately in 
front of the detector. X^ is assumed equal to n. 3. 

f The Gaussian local oscillator parameter Z is given by 
r ^ 

where w = 4- radius of L. 0. beam. It is assumed that 
w e 

equals 1. 

For a phase front tilt of X over the detector the heterodyne 
detector parameter, T, is found to be 'v 0.1. Note that, 

Kr 6 = K-^ = -|^ 4*’’^ "From Fig. 5. 3. 3(d) with Z = 1 and 
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FIGURE 5.3.3 


HETERODYNE EFFICIENCY CONTOURS 


Contours of equal valued heterodyne detection parameter -y as a function 
of Airy received signal parameter and Gaussian local oscillator para- 
meter 2^ ; (a) Kr^e = 0; (b) Kr^e = 1; (c) Kr^e = 2.1; (d) Kr^e = 3.1. 
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= 3 it cen be seen that Similarly, for a wave front 

tilt of the is found to be ~ 0.6. Using the information from 
Fig. 3.3 and the phase shifts calculated for the proposed ^-system 
in Tables IX and X the relative heterodyne detection parameter can 
be calculated. 

From Table IX it can be seen that the phase shift for the 

I I 

two extreme rays, Ray 1 and Ray 3 (for the existing /S -system) 
for returning beams on axis when the system is focused at 10 
meters is approximately 6.7 urn. It is also understood that this 
tilt takes on a shape that is given by Fig. 5. 3. 4(a) . This fourth 
power relationship permits a conservative approximation of the 
wavefront tilt at the secondary to be estimated at ~1.5 pm. 

This is further supported by the assumed Gaussian beam shape on 
the local oscillator beam and the Airy disk beam shape on the 
reflected beam. An observation of Fig. 5. 3. 3 indicates that the 
heterodyne detection parameter, Y , is approximately equal to 0.7. 

This performance tends to increase as the existing system 
is focused at longer ranges and approaches Y =0.75 for a focus 
at 30 meters. From Table X it can be seen that the phase shift 

I I 

for the two extreme rays, Ray 1 and Ray 3 (for the proposed 
/S -system) for returning beams off-axis when the system is focused 
at 10 meters is approximately 7.5 pm. Again the tilt takes on a 
shape that is given by Fig. 5.3.4 b). 

This forth power relationship permits a conservative approxi- 
mation of the wavefront tilt at the secondary to be estimated at 
2.0 pm. This estimate is also supported by the assumed beam shapes 
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of the local oscillator and the reflected/collected beari. An 
observation of Fig. 5. 3. 3 is that the heterodyne detection para- 
meter, > , is approximately equal to 0.65. 

This performance tends to increase as the proposed system is 
focused at longer ranges and also approaches the T = 0.75 for a 
focus at 30 meters. 
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4. surr,ARY 


Currently completed analyses of the two basic two-color ^-system designs 
which have been investigated, indicate that the two beam directing insertion 
mirror design will have the higher system performance efficiencyj consider- 
ing all other things equal. The study has resulted in a ray trace analysis 
which indicates that the existing system has a beam expanding telescope 
which couples with the Laser Doppler System to give relatively good heterodyne 
performance even when focused to 10 meters. 

The ray trace analysis indicates that the system can accommodate off- 
axis beam insertion without the serious penalty of substantially lower hetero- 
dyne efficiency. Relative performance between the on-axis beam and the off- 
axis beam may be undetectably small when measured in the field. 

The beam splitter combined two color /9 -system design has several desir- 
able attributes in that it requires only one detector, interferometer, and 
signal processor thus potentially becoming a smaller and lighter instrument 
package. This design may potentially be the more desirable, of the two in- 
vestigated designs, if performance efficiency for the system becomes less 
important and/or if simultaneous measurement of single particle signal events 
from Aj and A2 becomes more important. Further investigation and analysis 
of the data taken with the existing /3 -system will provide insight into the 
correct system design to choose for the two-color <5 -system. 
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7. Task 6: Evaluation of Applicability of Adaptive 

F i Iter ing to Pulsed LDV Signal Processing 

Adaptive filtering is appropriate for signals whose 
statistical properties are not known a priori, or whose 
properties change with time. In either case, adaptive 
processing constructs a Weiner filter matched to the 
particular signal. If the signal is stationary (not 
statistically changing with time), the Weiner filter is 
achieved and does not change with time. If the signal has 
non-stationary statistics, to be useful, the adaptive 
process must occur rapidly enough to "track” the changing 
signal properties. Since LDV signals are non-stationary 
narrow band random processes, it is appropriate to consider 
whether adaptive filtering can provide enhanced spectral 
estimation for LDV signals. While adaptive filtering 
applied to signal channel time series is reasonably well 
developed as "off-the-shelf” processing, applications to the 
spectral domain or to multi-channel data are not well 
developed. These latter applications would require 
significant research and development. 

1 ? 

Adaptive filtering has many forms, » each suited to 
specific signal processing problems. Forms originally 
considered for LDV applications were the Adaptive Noise 
Canceller (ANC) and the Adaptive Line Enhancer (ALE) 
implemented as an Adaptive Linear Predictive Filter (ALPF). 
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The ANC form is shown in Figure 6.1. This form requires two 
inputs: signal plus noise d(k) and a correlated noise 
source x(k) used as a reference channel. The reference 
channel noise is filtered and subtracted from the signal 
plus noise channel to improve the signal-to-noise ratio 
(SNR). Since a correlated noise source is not available for 
LDV processing, the ANC form is not appropriate. 

The ALE in the form of the ALPF is shown in Figure 6.2. 

There is a signal channel input to this filter, whose weight 
Wj are adjusted so that the e(k) output is a minimum (in the 
least mean square sense). The output x(k) is then the 
approximation to the signal without noise. This filter 
would be incorporated into an LDV system, as shown in Figure 
6.3, between the preamp and the final processor, if an 
additional processor is desired. In this way, the ALE is 
entirelt independent of the additional processor and may be 
removed, if desired, in the case of search and track 
failures. The purpose of the ALE is to improve the SNR into 
the processor, and it can be considered an add-on device. 
Uncertain at present are the convergence time (possibly too 
long) and tracking performance for a narrow band random 
process (aerosol signal). 
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Details on the ALPF: 


The filter coefficients w j , shown in Figure 6.2, which give 
minimum mean-square error for the e(k) are determined by 
values of the auto-correlation function of the input signal 
at various time lags."* These values of Wj constitute the 
•'Weiner filter". For a narrow band random process, this is 
a matched filter around the signal. As the signal changes, 
the filter must "adapt". To accomplish this, a "least mean- 
square" algorithm has been devised which continually seeks 
the "least mean-square" error."* This is a method of 
tracking the Weiner filter without calculating the auto- 
correlation lags. Important about this process is the time 
necessary for convergence for a narrow band random process, 
which has yet to be evaluated.^ 

Linear predictive filtering, itself, gives an estimate of 
the spectrum of the signal, called the "maximum entropy 
spectrum". Details of this are given elsewhere, ^ but, in 
summary, the spectrum is determined by the filter weight 
values. Figure 6.M shows the relation between the maximum 
entropy spectrum Q(w) and the transfer function of the 
"whitening filter H(w) which converts the signal x(k) in to 
white noise e(k). This method would be used only if the 
ALPF were the final stage of processing. 



Improvement in LDV Processing Offered by ALPF : 


Analytical formulae for evaluating the improvement in SNR 

offered by the ALPF in an LDV system are not available. 

However, evaluation of this effect can be accomplished by 

simulated narrow band random process data processed by an 

ALPF. The improvement on stationary data is shown in Figure 

6.5. This figure^ shows the velocity error to be expected 

in estimation of the mean frequency, as a function of SNR. 

The Cramer-Rao bound^ represents a theoretical upper limit on 

the possibility of improving the velocity error at a given 

SNR. The pulse pair (PP) curve is the theoretical response 

for the processing of 102^ points. This figure applies to a 

spectral width of 1 m/sec. Clearly about 12 dB of 

improvement above the PP curve is indicated. Thus it is 

reasonable to attempt to achieve this SNR improvement, 

although the processing techniques which would bridge this 

gap have not been identified. Also shown in this graph is 

the result for a linear predictive filter (LPF) as a final 

processor. A nominal improvement of about 3 dB is seen, but 

the performance is worse at the extremes of the curve. 

Improvement offered by a poly-pulse-pair processor as a 

final processor are not shown, but simulation has shown it 

to be very similar to the linear predictive filter curve. ^ 

Results for an ALPF followed by a pulse-pair processor or a 

poly-pulse-pair processor are not available. These results 

could be generated straightforwardlyt from existing 
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subroutines . 


Application of ALPF to Previously Recorded Data: 

An ALPF can be applied to previously recorded time series 
data in a straightforward manner. Power spectral data 
contain less information and are not amenable to 
conventional adaptive techniques. Attempts to apply 
adaptive processing schemes would require extending these to 
spectrum-based schemes, or to multichannel input adaptive 
filters. These are not ’’off-the-shelf” processing 
techniques, and would require development.^’^ 

Steps to improve mean frequency estimation based on power 
spectral data can be recommended. Noise spectrum averaging, 
from cases when no signal is present, provides a basis for 
appropriately weighting each channel. This correlation 
should be made. Then a circular mean estimation,^ indicated 
in Figure 6.6, should be applied as an inbiased estimator 
(in contrast to peak or linear mean frequency methods). In 
low SNR data, a plot of the circular mean frequency would 
dramatically indicate whether signal were present or not, 
since a noise frequency estimate would be found at any point 
around the circle, while a real signal would be consistently 
found in some angular interval. This would seem to be the 
best method of indicating the presence of real signal power 
spectrum data, since it is not biased by the noise spectrum. 
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Figure 6.2 

ADAPTIVE LINEAR PREDICTION FILTER 
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Figure 6.6 


CIRCULAR FIRST MOMENT FREQUENCY ESTIMATOR 




Summary : 


A potential gain of 10 to 12 dB in SNR beyond the response 
of pulse pair processing is theoretically indicated, but the 
processing techniques for realizing this gain have not been 
established. Therefore, it is reasonable to seek 
techniques, or combinations of techniques, which will 
approach this theoretical limit. Results have been 
presented which indicate that Adaptive Linear Predictive 
Filtering will contribute 3 dB of this gain, an improvement 
very similar to that offered by poly-pulse-pair processing. 
Combinations of processing techniques, such as ALPF and 
poly-pulse-pair processing, have not been evaluated. This 
could be accomplised with existing subroutines. 

Previously recorded power spectral data is not amenable to 
"off-the-shelf" ALPF processing, but would require some 
extension of existing techniques to the spectral domain, or 
to multichannel input data. However, noise threshold 
compensation measures and, in particular, the use of the 
circular first moment frequency estimates should give a good 
indication of the presence of signal. 
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8. Task 7: Correlation of Beta Experiment Data With LDV 


Data 


Software algorithms were developed to provide data 
processing capability to determine the backscatter 
coefficient from single particle measurements. These 
algorithms are capable of providing a correlation of beta 
measurement data with laser doppler data recorded in a 
different (volume) mode. Part of this work is also described 
in chapter 11. 

Single particle P predictions were made using data obtained 
an extreme focus of 40 meters. Figures 7.1, 7.2, 7.3, and 
7.4 show signal histograms of the data used to derive each 
prediction. The data used in the three figures was very 
sparse containing only 329, 705 and 192 single particle 
respectively. The correlation between predictions and this 
volume mode data (f = 40m) is not straightforward since 
these particles would be relatively large and the density of 
these large particles is also very low. 

Several assumptions were made in order to make these 
predictions. They are: 1) The noise level is constant at 

-59 dB. 2) The signal threshold was set at a S/N of 
approximately 7.5 dB. 3) The maximum signal bin is at -41 
dB and the correspondence between signal bin and signal is 

linear. Based on these values and other assumed system 
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parameters the beta inversion algorithms produced these 
results : 

1) Density of particles - 1 particle/m^ 

2) Average cross-section (backscatter ) - 33/‘^ 

3) Smallest cross-section seen - 1lM^(this is 
due to focus of ^0 meters) 


4) Beta = 2.2 

X 

C\J 

1 

o 

meter"^ 

for 

Figure 

7.1 

= 4.7 

X 

CM 

o 

meter"^ 

for 

Figure 

7.2 

= 1.2 

X 

CM 

1 

o 

raeter"^ 

for 

Figure 

7.3 


This system is currently being calibrated, however, these 
calculations are for an uncalibrated system. This fact plus 
the very small number of particles processed required that 
the predicted beta values be interpreted with caution. 

An additional case with a significant number of particles is 
contained in Figure 7.^. The time for this case is in error 
as shown, and the exact time is not available. Indications 
are that this case corresponds to a particle density of 
about 1000/m^. 
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The flight distances required to measure beta in the single 
particle mode are analyzed as follows: Since 


P 

and D 
where D 

'a 

A 

N 

inversion 

L 

we see 

L 

or log L 


= D CT / 4 ir 
= N/AL 

= total number of relevant particles/unit volume 
= average backscatter cross-section 
= LDV effective area 

= number of particles required for beta 
= flight distance/data set, 

= N a/Ur )3A 

= - log )3 + log (N ff /U r k). 


This shows that the flight distance is dependent on p as 
well as on the system and distribution dependent parameters 
A and a , and on the algorithm parameter N. Since A and F 
are averages over the backscatter probability distribution 
to be sampled, this distribution (or a set of expected 
distributions) must be known to predict L, independent of 
the prescribed p value. The accompanying Figure 7.5 shows L 
as a function of p for three different focal ranges, 10, 20, 
30m, and three different backscatter probability 
distributions: exponential, e"®*^ with B = .15 and power 

law, 1/<r® with B = 1.7 and 2,U, These distributions are not 
unreasonable, but do give a difference in flight distances 
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of almost an order of magnitude. Because they are different 

distributions, the densities for each case are different as 

A 2 

shown. Here, N was assumed to be 10 , the 1/e intensity 

radius was 5.5 cm and no beam truncation was assumed. 

Clearly, the smaller the parameter N^/A, the smaller the 
flight distance for a given /3 . For several reasons, the 
actual situation will be better than presented by these 
graphs : 

1. The 1/e radius is 4.6 cm giving a larger A. 

2. The actual beam is truncated, giving a larger A. 

3. N could be decreased as knowledge of the actual 
backscatter distribution increases. 

The actual backscatter distribution, itself, might give 
smaller ratios of F/A than represented here. 
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9. Task 8: 


Determination of 

The concept of arises when the return from a plane at 

the waist is used to infer the return from a volume 
distribution. A return from the waist plane can be written 

S„ = pj"r(x,y) dxdy 
waist 


where p is the reflectance in inverse steradians and F(x,y) 
represents the transverse beam profile. A volume return is 
written 


Sy = yS^F(x,y,z)dxdydz 


and used to define as 


Sy = PL 




gffji vA,y)dxdy 
waist 


Then if we ask what p value is implied by a measurement at 
the waist, we equate and Sy and find 

^ ^eff 


But we see that 
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Lgff = j F ( X ,y , z )dxdydz /^F(x,y)dxdy 

waist 

so that is really a geometric mean length determined 

only by the beam provile. 

An exact measurement of Lgff can be accomplished by making 
measurements in many planes and integrating. Let 

Sj„ = F(x,y)dxdy 
p plane p 

where closely spaced planes, separated by z, from signal 
beginning to signal end are chosen. Then, 

Sjn = .P S \ F(x,y)dxdy z 
p plane p 


.V _p F(x,y,z)dxdydz 

AZ 


So 


^eff - ^m 

and does not depend on the reflectivity of the target used 
for the measurement, or upon laser power. It does depend on 
correct optical alignment which is assumed. 
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1 . 


For the precise calibration of the system, the return from 
twenty-one planes separated by equal distances and centering 
on the waist should be measured. These planes should extend 
to points on either side of the waist such that the 
contribution from neglected planes determines the acceptable 
error. That is, if 1 percent error is acceptable, then 
measurement may be stopped when the next plane contributes 
less than 1 percent of the value of the integral Sjj^. The 
sum of these returns divided by the waist return times the 
plane' separation increment gives Lgff This set of 
measurements should be made at the operating focal range. A 
less precise value can be obtained by taking at least two 
sets of measurements at less than the operating focal range 
and extrapolating to the operating range. One can only 
guess at the precision of such a procedure since the real 
beam profile and its range variation are unknown. 

Experiments were made in January to determine Lgff for the 
NASA Beta System Laser Doppler Velocimeter (LDV). The 
effort was conducted using the system focused at 60 and 100 
meters. In December, measurements were made at these ranges 
but analyses indicated that additional data were needed. 
First the December measurements at 60 meters were made using 
a 10 inch aluminum wheel with 1-inch sulphur band on the 
periphery. These measurements were made out to distances 
approximately equal to two times the theoretical Lgj.^* . The 

p 

1/e beam size expanded to about one inch in diameter at 
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these distances and the beam consequently was larger than 
the sulphur band. Measurements were, therefore, made at 6Cm 
using the sand paper disk. 

Also in December, measurements were made using the vacuum 
chamber tunnel at 100 meters focus. Since the tunnel was 
only 100 meters in length, the laser beam was folded using 
mirrors to provide easy access to the beam at distances 
beyond the waist. Postion of the mirrors in the tunnel did 
not provide enough distance beyond the waist for a good 
estimate of Lgff Consequently, these measurements were 
repeated in January. 

Two characteristics of the spectrum analyzer were used to 
determine the sand paper wheel speed which was used during 
these measurements. A rotation speed sufficient to give 
Doppler shift greater than 300 KHz dead band, yet slow 
enough to keep the doppler band width well within the 
capability of the 300 KHz filter on the spectrum analyzer 
was essential. The shaded area in Figure 8.1 defines these 
limits. The Doppler bandwidth increases as the beam spread 
increases at larger distances from the waist. Consequently, 
to identify the presence of possible frequency truncation, 
measurements were compared using both the 100 and 300 KHz 
bandwidth filters. Differences between the signal levels 
measured in these bands were taken as indication of 
frequency truncation. A difference in signal level of two 
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dB was considered to be significant and data showing greater 
than this difference were not used in calculation of 
Slight differences in the signal levels were observed at 60 
meters indicating only a minor level of frequency 
truncation. However, significantly differences were 
observed in the signal levels at distances approximately 
equal to at 100 meters focus. Consequently, frequency 

truncation limited the distances to which measurements could 
be made at 100 m focal range. This resulted in an increased 
uncertainty in the calculated values of Lgff 

An additional comment should be made about the data 
analyses. Generally, signal levels taken at the end points 
exhibited about 18 dB attenuation from values measured at 
the waist. An exception to this occured at the 100 meter 
focus. For this focus, signal levels were attenuated only 9 
dB at the end points. The distance from the waist at which 
data could be taken was thus limited by the beam spread 
which resulted in a broader Doppler spectrum. The measured 
data was extended over a larger range for estimation of 

using a least square linear curve fit on the signal 
levels measured dB. This procedure resulted in a very good 
fit to the data as indicated by the variance and correlation 
coefficients obtained from analyses of the measured and 
curve fit data. The data was extended to approximately ^0 
dB attenuation levels using this procedure. Since the data 
experienced slightly faster than an exponential fall off. 
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this methodology results in an overestimation of . 

Results for the three ranges 30, 60, and 100 meters are 
depicted in Figure 8.2. These results are compared to 
theoretical values of Lgff untruncated obtained for a system 
with an aperture of .0^16 m. The majority of the difference 
between the curves shown in this figure can be explained in 
terms of over estimation of values by the procedure used to 
calculate and that the real beam is truncated. A curve 

for truncated ■ will occur between the two curves shown 

in the figure. As a result, the curves define boundaries 
for excursions within which the true values of l-gff 
truncated will fall. These cureves define the value of 
to about 20 percent in the worst case at 100 meters . 

In addition to developing the curve for Lg^.^. as a function 
of range, the beam profile was mapped for vertical and 
horizontal slices using 0.1 mm pin hole. Figures 8.3 and 
8.^ show the beam profiles at the waist for a 10 meter 
focus. These figures do not show any unusual features. 
However, the data given in Table 8.1 show evidence of 
astigmatism. 

Passive mapping of the beam profile was also made at 30 
meters focal distance similarly to the 10 m data. No 
unusual features other than the presence of astigmatism were 
found in this data. 
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Focus fmptcrs) 



Cross Beam Distance (Indies) 





Table 8.1 


Hcan- Wiv'tli fi'r Distantes fri>n; Waist at 10 Folus ir, 
Sht'wing the Effects of Astigmatism 


Posi tion* 

Vertical (inch) 

Horizonal (inch) 

Ih 

.2b'-) 

.2K7 

17 

.120 

.195 

18 

.255 

.150 

19 

.330 

.185 


* Waist = 16 feet on the optical bench. 
The heifjht of the beam is smallest at 17 
feet a foot earlier than the narrowest 
width of the beam at 18 feet. 
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10. Task 9: Calibration Using the Spinning Wire 


Data for the determination of beta in the single particle 
mode will consist of histograms of peak signal values. Such 
peak signal values will occur as the particle crosses the 
vertical central plane of the focal volume parallel to the 
beam. In calibrating the focal volume, sensitivity contours 
within this plane must be determined. These contours are 
approximately elliptical figures, each defined by a constant 
value of the backscatter signal to cross-section ratio, s/o . 
Therefore, mapping the intensity contours of this region by 
probing with a constant a probe provides a determination of 
the sensitivity contours. 

This mapping is accomplished with a rotating wire, and 
therefore, only specifies the s/o ratio up to a constant 
factor, since we do not claim to know o for the wire. (The 
wire is assumed to have the same "effective a " for each 
measured signal value). In order to understand the beam 
profile, contour plots can be made by taking data on 
transverse beam cross-sections. Then a subset of the same 
data can be used to obtain contour plots in the plane of 
interest . 
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Calibration Procedure: 


The following procedure is recommended as a first try at 
recording data adequate for contour plotting. It is very 
likely that attempts to process this data will require 
revisions in this procedure. 

1) The spinning wire should be configured so that the plane 
of rotation is parallel to the beam, intersecting the beam 
above the axis of rotation. 

2) The rotation plane should be moved through the beam in 
at least nine equal discrete steps with the fifth step being 
the beam misplane. Care should be taken to avoid backlash 
error. Steps 1 and 9 should be near signal-loss positions. 

3) Data should be recorded at each step as the wire passes 
through the vertical position, covering an angle of not more 
than about 30 degrees. This data should include the speed 
of rotation of the wire. 

4) A single test case should be obtained, consisting of 
many spectra from the same wire . position to investigate the 
effects of noise. If significant noise is found, several 
samples from each wire position will have to be averaged for 
the entire data set. 
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5) If possible, the same rotation speed of the wire should 
be used for each data set to avoid data processing problems 

6) A set of data files for each cross section should be 
recorded from beam waist to signal loss in each direction, 
utilizing at least 9 cross-sections as an initial trial. 

7) Recording of data in files indicating X and Z 
coordinates where X stands for wire position and Z for beam 
cross-section, from 1 to 9 would be very helpful in 
processing. 

8) Initial trial should be at 20 meters. Two subsequent 
calibrations should be made at other ranges to be 
determined . 

9) Software to process these data will have to register 
these data with respect to the beam axis, may have to 
average and interpolate data values, and will have to make 
contour plots. These plots will be used (integrated) to 
determine an area versus s/<r graph for input into the 
single particle processing algorithm. The location of this 
curve along the s/a will have to be determined by absolute 
calibration methods involving the single particle generator 
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Calibration Results: 


Ten sets of data were taken at a focus of 10 to 20 meters, 
respectively. The ten sets corresponded in both cases to 
from 1.0 meters in front of focus to 1.25 meters behind 
focus in steps of .25 meters. Backscatter signal was 
tabulated as a function of time as a wire (spinning on a 
wheel) penetrated the central focal area. Five time sets 
were taken at each position to minimize noise. The Fourier 
transform of each time set was calculated and then averaged 
to obtain one Fourier transform for each position in front 
of or behind the focus. Each Fourier transform then 
contains the information of backscattered signal versus 
height in the beam. This is true because each frequency of 
the Fourier transform corresponds to a distinct point in the 
laser beam focal volume (at this range). 

The method briefly described above is contained in 
Intergraph Corporation report 81-006 entitled, 

"Measurement System Design Study - Final Report (11 Jan. 

1982)«. 

Figures 9.1-9.10 show the Fourier transform at each 
longitudinal position for the 10 meter focus data and 
Figures 9.11-9.20 show the corresponding transforms for the 
20 meter focus data. In the figures shown, a 9 point data 
smoothing was employed to produce smoother plots. This data 
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smoothing did not appreciably change the features of the 
Fourier transforms, however, made data handling somewhat 
easier . 

Using 300 frequency bins about the center of each Fourier 
transform data set (i.e., each of the 10 range positions), 
contour plots of signal as a function of longitudinal 
position and vertical position (i.e., frequency bin) were 
calculated. Figure 9.21 is a linear contour plot of the 10 
meter focus data while Figure 9.22 is a logarithmic contour 
plot of the same data. Figures 9.23 and 9.2^ are contour 
plots of the 20 meter focus data. 

The 20 meter data is not ideal since the Fourier transform 
at some longitudinal positions is not complete. For this 
reason, the contour plots become unrealistic as the contour 
values decrease. 

Using these contour plots, area as a function of s/a can be 
tabulated. This is done by integrating the area within a 
contour of a particular value. Figure 9.25 shows a plot of 
area versus s/tr for the 10 meter and 20 meter focus data. 
The 20 meter data is more limited than the 10 meter data 
because of the previously mentioned problem. Also it should 
be noted that the data in both the 10 or 20 meter cases can 
be extended if data is taken from signal acquisition to 
signal loss along the longitudinal axis. The s/ a scale on 
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the horizontal axis is arbitrary as shown, the relative 
shape and spacing of the curves being the only important 
characteristics. Data from the single particle aerosol 
generator (see ARI Monthly Progress Report Number 9) is 
needed in order to establish absolutely the position on the 
s/cT axis where the "areas” begin. Noise is assumed to be a 
constant factor in Figure 9.25. 

Calibration Software: 

The software necessary for transferring these data to the 
Sigma 5 computer for analysis was developed. 

Two software routines were written for this project, and 
their documentation appear after the figures at the end of 
this task description. They are 

• A routine to collect wire wheel data and transfer to 
Sigma V. 

• A routine to simulate the remote terminal. 

The wire wheel routine reads data sampled by the Biomation 
Digitizer in sets of 20^18 data points from one pass of the 
wire through the beam. Up to 20 such sets can be collected 
before being sent to the Sigma V. The data is sent in ASCII 
format and as much may be placed in Sigma V file as desired. 
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This software has been developed to facilitate data 
processing . 

The remote terminal routine will allow communication to 
Sigma V from the test site. This routine will provide an 
easy means to verify the data transfer, process the data and 
determine the status of the Sigma in the case of transfer 
difficulty. Both routines have been successfully checked 
out using test cases. 

Additional software has been written to accomplish contour 
plotting of the focal volume profile from the wire wheel 
data. Following the previously mentioned documentation is a 
listing of a program which reads the data from a Sigma 5 
file and plot the contours. Following this listing is a 
contour plot of a test case sphere made with the contour 
plot subroutine, using data generated by a separate program. 
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I 


. SBTTL 

. MACRO 
. IF 

. ENDC 
. IF 

. IFF 


.ENDC 

.ENDM 

. MACRO 
. IF 


pACROG, 

ICALL 

MOV 

I 

NB,<. A> 
JSR 

.PRINT 

CALL 

PUSH 

NB,<:.A> 


.IFF 

.ENDC 

.ENDM 

. MACRO 
. IF 

.IFF 

.ENDC 


MOV 

.PRINT 

PUSH 

iPOP 

NB.C. A> 
MOV 

.PRINT 



. ENDM POP 

RETURN 

= RTS+7 

CR 

= 15 

LF 

= 12 

TAB 

= 11 

ESC 

= 33 

. Asect 

.LIST MED 

. ENABL AMA 
.ENABLE !lc 


' 1 , ; i ■■ 1 1 

I 

ASSEMBLY CONTROl. 

. A, .D 
'■.B",RO 

PC, ■ .f)' 

I 

I 

; MISSING ARG TO CALL MACRO 

.A 

•.A--,-(SP) 

; MISSING ARG TO PUSH MACRO 

.A 

(SP) + ,/. A' 

; MISSING ARG TO POP MACRO 


.•ESCAPE ASCII CODE 

; EXPAND MACROS 
; ABSOLUTE CODE 
; LOWER CASE 


1 




.SBTT 


IRTRP! HOLT 
RI1RP: HALT 


HrrRP! .WORD 


. WORD 


. WORD 


i WORD 


DUMY: 


1 


VECTORS 

0 


ILVEC 
IRTRP, 200 

RIVEC 

RITRP.200 

TRVEC 
BPTRP, 200 

PFVEC 

DPTRP,200 

600 

0 

2 

4 

# 200.6 
#2, 10 
# 200, 12 
#4. 14 
# 200, 16 
, 24 

#200,26 

1000 


;load interrupt section 



1 



1 


BTTL MAIN 
1 (jOO 

MAIN: MTPS I»200 :MI-‘.AI«Li: PROf. LSSi irv I N ILRi.iPTS 

MOV « DHLiR - 200 . . SP ; I M I 1 1 AL 1 Z E !: > T Ai: i 

MbV «20.,NDSPS ;S.IT « DATA SETC/SL I r f = 20 

MqV . «2048. , NSAPS P .IIT l» SAMPLEG/DA I A ’ .E r=^20Ai;:. 

< MOV «X3C1,R1 ;CHANNa. 1 XMIT TO TERMINAL 

CALL OUTSTR, ttCSOO •ANOUNlE PRESENCE 

CALL 0UTSTR,«CS02 TGIVE FILE NAME 

CALL 0UT3TR, ttCSOl 

CALL OUTSTR, ttCSOS 

CALL OUTSTR, IfCSO^ t PROMPT TO ENTER FILE It 

CALL NUMPR3 ; ENTER 3 DIGIT FILE' START It 

MqV R2iFILNUM ; SAVE FILE tt 

MOV ttXSCl.Rl 

CALL OUTSTR, «CS05 ? ENTER 2 DIGIT # DATA SETS/DEAM SLUI 

CALL NUMBR2 

MC|V R2,NDSPS 

INC NDSPS :ADD 1 TO ttSETS/SLICE TO MAKE COUNT lanUT 

* 

; ; ««« WAKE UP SIGMA SIGN ON «■ 

WAKEUP: MOV tt'H.Rl 

CALL STRSV,«CSA6 tWAKE UP SIGMA 

CALL CHECK ; CHECK FOR "H" OF HONEYWELL CP-V 

BNE WAKEUP ; TRY AGA IN IF NO GOOD 

MOV t»3.,R3 ;WAIT 3 SECONDS FOR COMPLETION OF SU-Mm I<l.SI urj:.! 

W3S mOv tfiooo. .r^ ;mili sec counter 

W2: MdV «238.,R". : •• TIMER 

wis sde R5.W1 

SOB R4 , W2 

SdB R3,W3 

SGNON: MOV «'*,R1 ; SIGMA RESPONDS WITH AT SION ON 

CALL STRSV,»CSA4 tSiON ON 

call check ; r; HECI^ FOR •• • RESPONSE 

BNE SGNON 5 TRY AGAIN IF NO GOOD 

; *♦****«# WIMI « « M II It II ItlHI H || H || n H |M| || mm 




MOV 

#XSC1,R1 

CALL 

OUTSTR, IICSOZ-. 

CALL 

CMETOF i 

STRT: MOV 

« 1 , DSCNTR 

MOV 

♦♦DDUE.DTAPTR 

MOV 

«DHDR,DHPTR 

MOV 

#XSC1,R1 

CALL 

OUrSTR, HCS07 

CALL 

NUMDR5 

MOV 

R2 . (<!DHPTR 

ADD 

tf2,DHPTR 

MOV 

HXSCl.Rl 

CALL 

OUTSTR, HCS08 

CALL 

NUMBRl 

MOV 

R2,eDHPTR 

ADD 

tt2,DHPTR 

MOV 

FILNUM, (5DHPTR 

ADD 

#2,DHPTR 

MOV 

#XSC1,R1 

CALL 

OUTSTR, «CS09 

MOV 

DHPTR , RO 

MOV 

«57.,R3 

MOV 

#RSC1,R4 i 

CALL 

INSTR 

DLOOP: MOV 

ttXSCl ,R1 

CALL 

OUTSTR, #CS 10 

DWAIT: CALL 

KB IN 

CMP 

R1 ,«'S 

BNE 

DWAIT 

CALL 

DREAD 

MOV 

ttXSCl ,R1 

CALL 

OUTSTR, ttCSl 1 

CALL 

0UTD3, DSCNTR 

MOV 

ttXSCl, R1 

CALL 

OUTSTR, «P I XL 

INC 

DSCNTR 

CMP 

DSCNTR, NDSPS 

DMI 

DLOOP 

MOV 

ttXSCl ,R1 

CALL 

OUTSTR, ttCS 12 

CALL 

DTRANS 

MOV 

ttXSCl ,R1 

CALL 

OUTSTR, ttCS16 

CALL 

KB IN 

CMP 

R1 , tf'N 


:i i., :i .. : :.i 

; SIGMA WAf .fr UP .?< SIGN ON COMPLi rF 
;COPY ME TO PILE 
; INITIALIZE DATA SET lOUNTER 
;SL1 DATA POINTER 
:Si:i DATA flEADER POINTER 

;ENILR •> DIGIT f.AMGPO M) 


; ENTER 1 DIGIT SLICE ft 


; ENTER COMMENTS (57 CHAR MAX) 
:&ET STRING DESTINATION 
?SET MAX It CHARS 
;SET STATUS REG CHAN 1 
;get comment string 

; READY TO TAKE DATA 
:GET CHAR ?< ECHO 

JWAIT IP NOT = "S" 

;get data set 

;SAY DONE 

IFOR DATA SET tt DSCNTR 

I 

t 

IDUMP DATA SET COUNTER 
5 CHECK rOR ff DATA SETS DONE 
ILOOP IF DSCNTR<--tNDSPS 

rSAY IiONE TOR ALL DATA SETS 
; TRANVT LR DATA 

IWANT TO CLOSE FILE 




rscF 


SGNOF 


JnIT 

^ E..:r> 

JMP 

STRT 

MOV 

ir!,Rl 

CALL 

STRSV, IH;:SA2 

CALL 

CHECK i 

DNE 

EscF : 

MOV 

#XSC1,R1 

CALL 

OUTSTR, «CS17 

CALL 

OUTSTR, ttCSOl 

CALL 

OUTD3.FILNUM 

MOV 

«XSC1 ,R1 

CALL 

OUTSTR, «PIXL 

CALL 

OUTSTR, HCS10 

INC 

FILNUM 

MOV 

ttXSCl ,R1 

CALL 

OUTSTR, ttCSl'I 

CALL 

KB IN 

CMP 

R1 , N 

BEQ 

SGNOF 

CALL 

CMETOF 

JMP 

STRT 

MOV 

#'C,R1 

CALL 

STRSV, 4»CSA5 

CALL 

CHECK 

DNE 

SGNOF 

JMP 

1 73000 


' ;Cu.... .V* r- n_il ir iJ:iT "iv^o ^ ^ 

; CONTINUE TO TAKE DATA OTHER WISE 

;SENU SIiiMA LGC F TO CLOSE FILL 

:1F SIOMA Lions NT GIVE " ! " TfcY AGAIN 

? SIGMA V 

; r I u: 

; « 


; CLOSED 

;dump fill: counter 

;D0 IT AGAIN? 


JCOPY ME TO FILE 


TSIGN OFF 


?GO TO KEYBOARD ODT (G?) 


1 










. SBTTL 


SUBROUTINES 


5 SUBROUT I NE STRS V 
5SEN|)S CHARACTER STRING TO SIOMA 
5 STRING ADDRESS SENT OVER IN RO 
;USE3 TRANSMITTE', SLRIAL CHANNEL 


V 


STRSV: 


1 $: 


BIT 

#200, XSC3 

BEQ 

STRSV 

TSTB 

< RO ) 

BEQ 


MOVB 

(R0)+, XDC 

BR 

STRSV 

RETUF 

N 


;TEST DIT 7 IF SET XMITTLR RI.ADY f tJI XT 
PIT 7 CLEAR? WAIT 
TEST STRIN'i CHAR 
EXIT IF NULL 
SEND OUI CHAR 


I IIG 
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1 


1 


1 


U lEC K i 
3t: 


2 *! 


1 *: 


5«: 


SLiniROUTINL' CHLCK 

CdhtPAREo CHAR IN fil TO lMAI< ON CERIAL INPIJI 
Tll]j oEE IP c.ICMA V OlOf C C.i.iFiREC I R’LCPONSE T'j C OliMANP'' : 

CI|Ar< TO E^E COMF AInLLi CENT IN R1 

iJ Rl.TSTCHR :SAVr; TEST EUAEf 

; II •cLCONliS TO WAIT mi.; IvrS.^PONSE 
; MILL ISEC rijI.ilMTl R 
TLOAM J MS I I MLR 

;'.n[ :i;i;: in t 7 i.i ian ivt cv. 

; IE SET CIIIXI: EIIAR 

JWAIT 1 MS THEN Rlj: HE.i: 1^ |\i:CV. DONE 
TCi.'UMT J OT'O ms E0T\ 1 SEC 
: COUNT II SECS 

;PRNT •■’SIOMA POLS NOT RECr ONP" 

;SLT Rl^l SINCE NO GOOD 


MO 
MOC| 

Mod 

MoO 

bit! 

PNG 

SOB 

SOI^ 

SOB 

CAUL 

Mov! 

PR I 
MOVB 
Bid 
CMP 

RETlURN 


R1 , TSTCHR 
ttlO. ,R.S 
ttlOOO. . R-^ 
«238. .RS 
ttSOO, RSC3 
It 

R?,2t 
, 3t 
R3, 

STRNLP, IICSAl 

ttl,Rl 

5t 

RPC3, RO 
«1 77600, RO 
RO , TSTCHR 
3t 
R1 


; GET Cl lAR 

; CHECK FOR CORRECT CHA|-C 
TIP NOT CONTINUE WAIT 
TIP GOOD SET Rl^O 


T SUBROUTINE STRLP 

’S^NDS CHARACTER ;5TRIN0 TO L.INE PRINTER 
TStRINO ADDRESS SENT OVER IN RO 
;U3ES TRANSMITTER SERIAL CHANNEL 1 


TSAME LOGIC AS STRSO 


STRNLPS Bn 


#200, XSCl 

DEC 


STRNLP 

tsi;b 

( RO ) 

BEO 

It 

MOV 

D 

(FVO) + , XPCl 

PR 


STRNLP 

Its return 



1 ' 1 


Ell iHI. 



] 


1 




K3 

(jJ 


nUTGTR: 

21:: 


1*: 


INSTR: 


!♦: 


2ti 




! " :i, 1 ]_ i i i 

I 

; SUPROUt|i NE OUTSTR 

{OUTPUTS! STRING ON SERTOL LINE 

{WHOSE TRANSMITTER STATUS API it<E: S'.. IS IN R1 

{STRING AULiRr:SS IN RO 


{ USES 

MOV 

ADD 

DIT 

DEQ 

TSTB 

DEQ 

MOVB 

BR 


fvOi, R1 , R2 
Rl ,R2 
H2. ,R2 
H200, (RJ ) 
2t 

(RO) 

It 

j( RO ) + , ( R2 ) 

Et 


RETURN 


{PUT AMIHVCSS Ml' XMIT DATA DMI f I R IN Ia 

{TE-:.T nn / - IE SI. I XMIIH.R l\T (UiY En 

{Dir 7 TErOR WA(T 

{TEST :,:.TRINi:i Eufc NULL 

{[.XIT IE MULL 

{SEND CIIAIvACTEIV 

{ Loot'- 


{ SUBROUTjl NE I NSTR 

{INPUTS BTRING ON SERIAL LINE 

{WHOSE RfeCIEVER STATUS ADDRESS is IN R4 

{STRING DESTINATION ADDRESS IN RO 

{MAX # CHARS IN R3 

{USES R0Lr2,R3,R4 


{ CALLS KD I N , OUTSTR 
MOV R4 , R2 

ADD |»2,R2 

MOV ho,PTR7 

MOV R3.PTR6 

CALL kBIN 

MOVB Rl,(RO)+ 

DEC R3 

TST R3 


{ADDRESS or RECIEVER DATA 
{SAVE Fa(.» in case comments 
{SAVE MAX IT CHARS 
{ GET .!< ECHO CHAR 
{ SAVE 

{COUNT CHAITS INPUT 
{IE “ MAX It TRY AGAIN 


DEQ 

CMP ^i,«o:r> 

BNE it 

MOVB feo,(RO)+ 

BR i5t 

MOV HXSCl.Rl 


{CHECK EUR CARfVAGL RETURN 
{PUT NULL AT END 


CALL 

MOV 

MOV 

BR 

RETURN 


pUTSTR, ttCSlS 

PTR7,R0 

N‘TR6,R3 


COMMENTS TOO LONf. 
RESTOr^lE RO 
" R3 

GET COMMLNTS AGAIN 


liiJEEER IN R 
TOO LONG 


] 1 ■ ■■ 1 


HI >: I Dr II. 



c hinror: 


DTRANS: 


CMTS! 


TLOOP: 



1 1 1 

x'NE .-.l•.£TOr ' 

1 

1 ) "1 

1 ■ 

; SUNOS 'i 

a CM A "COPY ML TC 

FILL 

F ILC NAME II Hit" 


;USE3 rJ 





; CALLS ■: 

■ilRSV , niJTDS , CARRET 



MOV 

«'.,R1 




CALL 

STRSV,«CSAO 

; COMMAND SIi iMA "COF ’Y Ml' fO 

f I LL " 

CALL 

STRSV, ttCSOl 

; 3I:ND 

FILE NAME 


CALL 

0UTD3.F-ILNUM 




CALL 

STRSV, «PIXL 

; SI' NO 

FILr NUMIfER 


CALL 

CARRET 

-, St ND 

<CR> ?<■ CHECK FOR 

REGION 

RETURN 





? SUBROU" 

INE DTRANS 




; TRANSFE 

:RS DATA TO SI DMA 

V 



? CALLS 

^TRSV,0UTD3,0UTD1 

, OUTD!; 

, CARRET , CHECK , L 1ST , OI IT3TR 

? 0UTD2 




MOV 

«DHDR,DHPTR 

; RESET 

DATA HEADER POINTER 


CALL 

OUTDS.GDHPTR ! 




CALL 

STRSV, #P I XL 

? SEND 

RANGE 


ADD 

»2. , DHPTR 

; BUMP 

HEADER POINTER 


CALL 

OUTDl,eDHPTR 

» 



CALL 

STRSV, «PIXL 

; SEND 

SLICE II 


ADD 

#2. , DHPTR 




CALL 

0LITD3,GDHPTR 




CALL 

STRSV, #P I XL 

; SEND 

FILE « 


ADD 

#2. , DHPTR 




CALL 

STRSV, DHPTR 

; SEND 

COMMENTS 


MOV 

«•'. ,R1 




CALL 

CHECK 

; CHECf 

FOR " . " FROM S I GMA 

V 

BNE 

CMTS 

;TRY AGAIN IF NO GOOD 


MOV 

#1,DSCNTR 

;SET DATA SET COUNTER=^l 


MOV 

«DBUF, DTAPTR 

;SET DATA POINTER 


CALL 

0UTD2,D3CNTR 




CALL 

STRSV, #P I XL 

; SEND 

SETIt 


CALL 

CARRET 

*, SEND 

<CR> S< CHECK FOR " . " 


MOV 

DTAPTR.PTRl 

• 

» 



MOV 

NSAPS.PTR3 




CALL 

LIST 

; SEND 

DATA SET 


MOV 

PTRl .DTAPTR 




CALL 

OUTDS.CDTAPTR 




CALL 

STRSV.ttPIXL 

; SEND 

WIRE GPETD 


CALL 

CARRET 


;sEND ;cr;: check 

FOR 

ADD 

#2,DTAPTR 

; BUMP 

DATA POINTER 


CALL 

OUTD2, DSCNTR 




MOV 

ttXSCl ,R .l 




CALL 

OUTSTR, 1fC313 


• traM' I'Te; pour n iK' 

DA 1 A • 



' 1 


1 




CALL” ' 
INC 
CMP 
DM I 

RETURN 



. BLKW 

? READS 
; PACKS 
!USES F 

,0 


dread: 

MOV 

h- * 

N5 ^ 


MOV 

Ln 


1*; 

TST 

'• 


BPL 


2t: 

BIS 



BIC 


3t: 

TST 

BPL 

MOVB 

DEC 

BNE 

MOV 

MOV 

SWAB 

BIC 

MOV 



ADD ! 

RETURN i 


OUTSTFv',«PlXL 

DSCNTR 

DSCNTR.NDSPS 

TLOOF- 


; II 

si:r LOUNniR 
iLOdP IP Ii:.i;.NiR' ;Nlr..l 


1 . 


BIOMAT I ON DIOITIZPFn- EiATA (; 
2 BYTES INTO 16 EFIT WORP 
0,R1 


I'lT BYTF.S) 


nsaps.ro 

DTAPTR.Rl 

DRCSR 

1 $ 

# 1 , DRCSR 
#1, DRCSR 
DRCSR 

3* 

DRIN, (R1 )+ 

RO 

It 

R1 . DTAPTR 
DRIN, WSPEED 
WSPEED 

#177400, WSPEED 
WSPEED, 6DTAPTR 
#2, DTAPTR 


:SET BYTE COUNT 

; storaof: array 

;TEST REOUEST D 
: BRANCH IF SIGN BIT CLP Ale 

TOGCiLE 'CMD' LINE 
LINE 

TEST roi^ REQUEST B 
HIGH AGAIN 

: STORE DATA 
? DECREMENT BYTE COUNT 
:1F NOT ZERO, DATA REMAlNl 
:SAVE R1 IN DATA POINTER 
IREAP WIRE SPIED IN MOST 5 
VSWAI-’ BY FES 

;AND CLEAR DUS 8 - IS 
;SAVE WIRE SPEED 
;DUMP DATA POINTER 


SIGNIFICANT BYTE 



L 1ST: 
1 ♦ : 




4*s 


N3 


? 7-’DrOL ' 1 i ^ 1 

; SENDS hIaTA to sigma V IN ASCII 


; UNPACKS 

DATA PUEFLR 

FOR 0 HIT VALUES TO SEND 

IN S DIGIT A 

CLR 

,PTR2 

;[£NlLf^ WITH PTR1=:DA1A 

AtiDRESS 

MOV 

|#15. ,PTR4 

;AUP PTR-;:.:r^« ITEMS 


CALL 

0UTD4.PTR2 

; PRINT TIEM NUMBER 


CALL i 

iSTRSV,#PIXL 

ISLND DATA ITEM II 


CLR 

BYTEF 

; CLEAR BYT E FUYT 


INC 

;PTR2 



INC 

DYTEF 

ISET - 1 SINCE UUlPLiriNG DYIE 

Move 

ePTRl .BYTE 



CALL 

0UTD3.BYTE 



CALL 

STRSV, ttPIXL 

JSEND DAI A riLM 


CLR 

BYTEF 

: CLEAR BYTE FLAG 


ADD 

«1 .PTRl 



CMP 

PTR2.PTR3 

;enuf items? 


DPL 

4t 



DEC 

PTR4 



DPL 

2 % 



CALL 

CARRET 

JSEND CARRAGE RETURN 

CHECK FOR " . ’• 

BR 

It 



CALL 

CARRET 



RETURN 




; SUBROUl 

INE CARRET 




?NO ARGSt 
CARRET! MOV 

calL 

CALL 
BNE 

RETURN 


T SENDS CARRAGE RETURN TO SIGMA CHECKS TOR RESPONSE 


PASSED 

,R1 

STRSV,«CSA7 

CHECK 

CARRET 


;SEND CR 
; CHECK FOR •• . *• 

;TRY AGAIN IE NO GOOD 


BLKW 


10 
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OUTDl! 


OUTD2: 


0UTD3: 


OUTD4S 


OUT DSP: 


0LITD5: 


DIGIT: 


WORD: 


. SDTTL 


terminal I/O ROUT I Nr:; 

N-DIOIT LIECIMAL UUIIUI KOUTINL:. FUlc 1 N 1 1. lx'l<UI T U 


MOV 

MOV 

BR 


« 1 , RS 

ma.R4 

DIGIT 


; 1 - D 1 G 1 T LILC 1 M AL OU 1 1 II T 


MOV 
MOV - 
BR 

MOV 

MOV 

BR 


«2,R5 

«K10,R4 

DIGIT 

#3,R5 

«K100,R4 

DIGIT 


T2-IUGIT DECIMAL OUTI'IJT 


; 3-DIGIT DECIMAL OUTPUT 


MOV 

MOV 

BR 

MOV 

MOV 

BR 


«4,R5 

«K1K.R4 

DIGIT 

45, R5 

#K10K,R4 

DIGITP 


MOV #5.R5 

MOV maoK,R4 

JSUBRpUTINE DIGIT 
ITEST 
MOV 

TST BYTEF 

BEQ WORD 

TSTB RO 

BPL DIGITP 

NEGB RO 

MOVD (R3)+ 

BR DIGITP 

TST RO 

BPL DIGITP 

NEG RO 

MOVE #-'-,<R3) + 


: 4-DIGIT DECIMAL OUTPUT 


; 5-DIGIT POSITIVE OUTFUT 


: 5-DIGIT DECIMAL OUTPUT 


BYTEF FOR BYTE OR WORD PUTS IN PI XL IF NEG 

#PIXL,R3 

;TEST BYTE FLAG 1=BYTE 0=WORD 
;WORD IF 0 

;test byte 

? BRANCH IF + 

;PUT IN PI XL 

;TES1 word 
T BRANCH IF + 


; PUT 


IN PI XL 
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MIIV I lK‘f ),+ ,,Ki ••'Al.l !•»/ li fn .1 N I * • ' , 

2f. ' -^r.r' ■ .k),Kv fill 1,1 I 1.1 1 1 

BLO 3 % 

INC R2 

SUB Rl.RO 

BR 2t 

3t: ADB «"0,R2 

MOVB R2.(R3)+ 

SOB R5.DIGITP 

MOVB It ' , ( R3 ) 4 

CLRB ( R3 ) + 

RETURN 

PIXL: .BLKW 3 

0000 

KlOK; 10000. 

KIK: 1000. 

KlOO: 100. 

KIO: lO. 

Kl: 1. 

? N-DIGIT DECIMAL INPUT FOR INTEFIRUPT HANDLERS 

NUMBR 1 : MOV t» 1 , R3 

BR I NUMBR 

NUMBR2: MOV «2,R3 

BR I NUMBR 

NUMBR3: MOV «3»R3 

BR I NUMBR 

NUMBR4: MOV t»4.R3 ^ 

BR I NUMBR 

NUMBR5! MOV #5,R3 

INUMBR: CLR R2 ; ENTER W/R3^==ll DIGITS 

1«: MOV R2.R1 

MUL <»10.,R1 

MOV R1.R2 

CALL NMBR 1 ; GET A DIGIT 

ADD R1.R2 ;AND ADD IT. 

SOB R3,l* 

RETURN ; RETURN W/tt IN R2 


BLKW 


1 
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.1 


I 


. SDTTLi UTILITY INPUT r\nuriNL:.. 


CHOR: 


CALL 

RETURN 


KB IN 


;OLT A CUAf< FROM USER 


NMBRi: call 

KB IN 

ADD 

#-/-0,Rl 

BMI 

NMDR 1 

CMP 

#11, R1 

BLO 

NMBRI 

RETURN 



;OLT A ItlOlT I ROM USER 

;LLS-3 TUAN ZERO 

:I3 IT REALLY A DIOITT’ 

;no u.nort: 


LETR: 


CALL 

CMP 

BMI 

CMP 

BPL 

. RETUR^ 


KB IN 


#172, R1 
LETR 
#100, R1 
LETR 


;GET A CHAR A-l' ONLY 


; RETURN VALUE IN R1 


GET A CHAR FROM USER (KB INTERRUPT DISABLED) 


KB IN: 

MOVB 

TKB,R1 

; EMPTY BUFFER 

2*i 

TSTB 

TKS 

;WAIT 

FOR FLAG 


BPL 

2*. 



KBINl! 

MOVB 

TKB,R1 

;CHAR 

NOW IN R1 


BIC 

#177600, R1 



Iti 

TSTB 

TPS 

; ECHO 

IT 


BPL 

1« 




MOVB 

Rl.TPB 




RETURF 





.BLKW 

1. 





SBTTL 


CGUNTER3 , POINTEK'S VAKlAliLEIS 


i COUNTERS 


ND3PG: 0 

NSor-3; 0 
riLNUM: 0 
DSCNIP; 0 


. BLKW 


20 . 


; POINTERS 


DATA SETS/GLICL IN DEAN 
; ttSAMPLES/DATA GET 
;F1LE U Eul\- SIGMA V (S DlulTS) 
5 DATA SET COUNTER- 


(jj 

o 


DHPTR: 

DTAPTR: 

PTRl: 

PTR2; 

I-TR3: 

PTR4: 

PTR5: 

PTR6: 

PTR7: 


0 

0 

0 

0 

0 

O 

0 

0 

0 


JVARIABL 


PS 


BYTE: 0 

TSTCHR: 0 
BYTEF: 0 ' 

WSPEED; 0 


;data header pointer 
;data buffer pointer 

; GENERAL-PURPOSE POINTERS FOR LIST 


J TEMPORARY STORAGE EOR L r3T 
5 TEST CHARACTER* FOR *SUB. CHECK 
»BYTE FLAG FOR BINARY TO A*:;CII RAIITINE 
JWIRE SPE.ED 8 BITS 



r 


u> 


CSOO: 
CSO 1 : 
CS02: 
C?03: 
CSO^ : 
C:'0?: 
CS06: 
CS07: 
CSOS: 
CS09S 
CSlOs 
C31 IS 
CS12: 
CS13S 
CS14: 
CS15: 
CS16S 
CS17S 
CS18S 

CSAO: 
CSFILE! 
CSAl 

csa: 

CSA3: 
CSA4! 
CSA5! 
CSA6! 
CSA7! 
CSA9! 
CSB2! 
CSD3: 
CSB5! 
CSB6! 
CSB7! 


. ASCIZ 
.ASCI 2 
.ASCIZ 
.ASCIZ 
.ASCIZ 
.ASCIZ 
‘.ASCIZ 
.ASCIZ 
.ASCIZ 
•ASCIZ 
.ASCIZ 
•ASCIZ 
•ASCIZ 
.ASCIZ 
.ASCIZ 
.ASCIZ 
.ASCIZ 
.ASCIZ 
.ASCIZ 
•ASCIZ 
•ASCIZ 
•ASCIZ 
•ASCIZ 
•ASCIZ 
•ASCIZ 
•ASCIZ 
•ASCIZ 
•ASCIZ 
•ASCIZ 
•ASCIZ 
•ASCIZ 
•ASCIZ 
•ASCIZ 
•ASCIZ 


I liiANSI'LR flvUuRAM RIInliY" 


SIOMA V" 


^CR><LF>"WIRE WHQIL DATA COLl.ti f inN ANIi 
•fW I RE WHEEL" 

CRXLF>"FILE NAMES ON STOMA WILL PE"<CR;: 'ILF > 

FOLLOWED BY 3 DIGIT DECIMAL IT" 

CRXLF>"PLEASE ENTER S DIGIT FILE START TT" CIE CLP >"« 
<iCRXLF>"EN rt^fi 2 DIGIT IT OF DATA SE T S/DLAM SLirr - " 
<CRXLF>"SIGMA. V AVAILABLE FOR DATA TRANSFER" 

•CCRXLF:>" ENTER 5 DIGIT 02767 MAX) RANFF O H) — " 
<CRXLF>"ENTERi I DIGIT SLICE « " 

CR><LF>''ENTEr' FILE COMMENTS (57 CHARS MAX ) "0;R><Lr> 
CRXLF>"READY TO TAKE DATA ENTET^ - S - " 

CRX;LF>"FINISHED taking data for DATA SET tt " 

CRXLF>"FINISHrD FOR ALL IjATA SETS - TRANSI LIVING 10 
CCRXLF>"FINISHED TRANSMISSION FOR DATA SET « " 

:CRXLF>"DO IT AGAIN (Y)ES OR (N)0 " 

CCRXLF>" COMMENTS TOO LONG TRY AGA IN"<CR>:LF>" w " 

:CRXLF>"WANT TO CLOSE SIGMA V FILE (Y)FS OR (N)nr^ •• 

:CRXLF>" SIGMA V FILE " 

' CLOSED ! " 

"C ME TO '• 

'SIGMA FILE NAME 

:CRXLF>"SIGMA V DOES NOT RESPAND" 

:esc>"F" 

:CRXLF> 

' B I LBRO . EC32505 " <CR> 

•OFF"<C;r> 

CESCXESO TEQUIVALENTE TO BRAICE 

:CR> 

fcCR><LF>" ENTER SIGMA V FILE NAME (11 CHARS MAX) THEN CR" 0';R; T ( 
fCCRXLF>" USING SIGMA COMMAND" :CRXLF> 

jCCRXLF>"FILE NAME TOO LONG TRY AGAIN PLEASE" CCRXLF:; 
|CCRXLF>"WANT TO RUN AGAIN'' Y(E'S) OR N(0)" 

|CCR::><LF:>"FILE name ok? Y(ES) or N(0)"CrR- <| r-: 
ccr:.><:lf>"enti£r aoain"<cr;: lf> 



SBTTL storage; AL LOCAT I ON 


« 3 - 

+ 

3 

t- 

(/) 

C. 




> 

5 


2 

.111 



3 3 

- - V*- V 


o 

XU 

> 

LU 

• 

II 

II 

* A 

a a 

• • 

II 

2 

U) 

• 


3 


•• «a 

a 



£C U. 

2 

CO 


Cl 2' 

U 

c 


X a 

»- 


a 

a o 

Cl 
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1 ’) 


TITLE DUMBT V 1.0 1/21 /.I:; 


. SBTTL I DEN I I E AC T I ON 


» *« W *MM» » » « tH! tl « im M « * mm (HH; H « « II l( tf 1» W « If B 

DUMBT: DUMB TERMINAL ROUTINE TO TALK TO SIOMA V 

TRANSFERS 10 SIOMA V 


WRITEN FOR NASA/MARSHALL 
BY IVAN BURROUOHS/APPLIED RESEARCH INC 

131 LONGWOOD AV. 

HUNTSV I LL E , AL SSSO'l 
(20S) 033- 6907 

; -IHHHHHfi-IHflMfWfffHffHf lf*B tt#mfHf mf #■« ftB If B If «ff Ifff ff fl ffBff B » II H ft If 




RO 

R1 

R2 

R3 

R4 

R5 

SP 

R6 

PC 


.SDTtL niRFINITIONS 
; PROCESSOR IJEr 1 IM I T I ONS 


es 


tr 


•/.o 

y.i 

y .2 

y.3 

y.A 

y.Tj, 

*/.6 

3P 

7.7 



SERIAL LINE UNIT 


oj 


TKS 

TKB 

TPS 

TPB 

KBVEC 

PRVEC 

RSCl 

RDCl 

XSCl 

XDCl 

RSC2 

RDC2 

XSC2 

XDC2 

RSC;3‘ 

RDC3 

XSC3 

XDC3 

R2VEC 

X2VEC 

R3VEC 

X3VEC 

DRCSR 

DROUT 

DRIN 


177560 
TKS+2 
TKS+4 
TPS+2 
60 
64 
TKS 
TKD 
TPS 
TPB 
1 76500 
RSC2+2 
1 76504 
XSC2+2 
176510 
RSC3+2 
RSC3+4 
RSC3+6 
300 
304 
310 
314 

167770 

DRCSR+2 

DROUT+2 


;KB STATUS USES CRT KEYBOARD 
!KB BUFFER 

;PR STATUS CHANNEL 1 SERIAL XMIT TO CRT 

;PR BUFFER 

;KI) INTERRUPT 

;CRT INTERRUPT CHAN I 

tRECIEVER STATUS SERIAL CHANNEL 1 CRT CI YDi lAMi 
? *' DATA BU FER " " 

; TRANSMITTER STATUS SERIAL CHAN 1 " 

; ’■ DATA BUFFER 

5RECIEVER STATUS SERIAL CHANNEL 2 LP KEYBi iC.Rti 
i '• BATA BUFFER " " 

; TRANSMITTER STATUS SERIAL CHAN 2 FOR LINT PICINTPr< 
? " DATA BUFFER CHAN 2 

iRECtEVER STATUS SERIAL CHAN 3 FOR SIGMA V 
? •' DATA DUFFER’ " " 

: IRAN'SMITTER STATUS CHAN 3 
' " DATA DUFFER CHAN 3 

TRECieVER CHAN 2 INTERRUPT VECTOR LP KM 
IXMIT IT-'R CHAN 2 INTERIiUP T VLCiTOR 
•.RECILVER •’ 3 

TXMITIER " 

JDRVll COMMAND STATUS REGISTER 
; OUTPUT DUPFEIi- 
; INPUl DUFPER 


This Page Intentionally Left Blank 



SBTTL 


IFxTRP: 
HI Tlvp: 

BPTRP: 


DLU1V: 


MAIN: 

Luup: 

L-.. I 


VECTORS 


O 


HALT 

HALT 


- 

ILVEC^ 

.WORD 

I RTRP . 200 

:s 

RIVEC 

. WORD 

R I Tf<P . 200 


TRVEC 

. WORD 

BPTRP , 200 


PFVEC 

WORD 

DPTRP.200 


600 

CLR 

0 

CLR 


CLR 

4 

MOV 

«200, 6 

MOV 

«2. 10 

MOV 

«200. 12 

MOV 

H4> 14 

MOV 

»200, 16 

MOV 

114 . 24 

MOV 

, 41200,26 

DR i 

1000 


I 


;LOf^D INTERRUPT SECTION 


I 


.SEHTL 
; S I GMA 
; CRT 

MTPS 

MOV 

CALL 

CALL 

CALL 


rU'VfN 

V ON SEIvlAL CHANNEL 0 


1 000 
! 11200 
I HXSCl.Rl 

i OUTSTR. HCSOCi 
I KB IN 
■ SVIN 

i i. 


1 

; DISABl i: l -Ri ii rSSOR IN f ERUI-'TS 
; CHANNEL 1 XMir TO TERMINAL 
; ANCiUNl E E-RESENCE 

iCET CHAR. SENG TO SIGMA S 
;GET SIGMA CHAR PRINT 

I . .. I -..; L -_ f _ 


ECHO 

I 
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•')SUrf^'''ViT 1“' 
?GET3 


SVIN: 




■iSVT 

tHAR KROM 


1 


1 *! 


TSTB 
BPL 
MOVB 
BIG 
TSTB 
BPL 
MOVB 
RETUrtN 


:;I0MA 

RGC3 

It. 

RBC3.R1 
«1 77600. R1 
TPS 
2t 

Rl.TFB i 


~\ -~] '^'1 1 

IF PRESENT S< PUTS TO CRT 

;test if char present 

;EXIT IF NOT 

;CHAR IN R1 
;TEST CRT STATUS 

;SENi:i TO CRT 


1 1 


{SUBROUTINE KB IN 

?GET^ A CHAR FROM CRT IF PRESI'.NT PUTS ON CtcT 


KBIN: 
KB I Nil 

ntr 

It: 

EX! 


TSTB 
BPL 
MOVB 
BIC 
TSTB 
BPL 
MOVB 
TSTB 
BPL : 

, MOVB 
RETLH|?N 


TKS 

EX 

TKB.Rl 

« 1 77600, R1 

XSC3 

4t 

R1 , XBC3 

TPS 

It 

R1,TPB 


? IF NO CHAR IN DUFFER EXIT 

;CHAR- NOW IN R1 

{TEST STATUS ON SERIAL CHAN 3 
{WAIT IF NOT F^EADY 
{SEND TO SIGMA V 
{ECHO IT 


{OUTPUTS STRING ON SERIAL LINE 
{WHOSE TRANSMITTER STATUS ADDRESS IS IN R1 
{STR:NG ADDRESS IN RO 
{USES R0.R1.R2 


OUTSTR! 

2t! 


It! 


MOV 
ADD 
BIT 
BEQ 
TSTB I 
BEQ 
MOVB I 
BR 

return 


R1.R2 
tt2. .R2 
«200.(R1) 
2t 
(RO) 

It 

(R0)+. <R2) 
2t 


{F-UT ADDRESS OF XMIT DATA BUFFER IN R2 

{TEST BIT 7 - IF SET XMITTER READY FOR Nl.XT FiYir 

{BIT 7 clear wait 

{TEST STRING FOR NULL 

{EXIT IF NULL 

{ SEND CHAFvACTER 

{ LOOF'' 


CSOO! .ASCkz <cr><;lf>"dump terminal FVEADY ! ! " <cr> :lf> 

.END' DMMY 



Contour Routine and Sample 



0'5:3-? l£C 23/'61 DC/?K jC 6 , t: I 


C L' ' Vf J T “ C 2 0 -J 2 J .7 T.I C 2 0 d J 

C0'^-'C‘-/Si'?/®(1012) 

LI ‘£‘.S10\ X(9), Y(10n,S(9,101) 

C-LL r^ITT(9o:0) 

FL^ro 5 
D^ = .'.953 

i.r=i^i 

\5 = o 
JE s K S 

•.n = 2C24 

iN r P s 2 k- 5 
Mj=l 1 

FSs 20.£tQfa 

Trt5.0035 

f,D2 = ND/2. 

DO 7 L=1/NDP 

TRCL3=0. 

7 TI(l)=0. 

LS = :iS/2iJ 

DO 1 1=1, LS 

1_ RE A DX5_f JJ) 0 0 J _L T P C L J.7 L =1 , 2 0 2 9 ) 

CfiLL FFT(\D°,-vJ) 

DO g L=1,ND 2 

4J P(L) = (TR(U*T«(L)+TI (L)*TI (L) )/(FS*FS) 

W^ITE(fe,50Q0) C P(L)7L = 1,1012) 

1 = 0 

CALL. SERCh(\wi,N02,I,PT.) ) 

•HKlTEClOI ,20001 NMl 
0OTPuT..P_TI 

LL=1-ND 


PEkWND 5 

_.DO_.fe_I=J/fMS 

DO 9 L=l, NDP 

TPCLJiO. 

9 TI(L)=0. 

PUT 


READ(5/1000) (TP(L)7L=1,202«) 

tt PII £ L^.JLQS> 0) (TRt I ) y L s 1 7 2.0 2JU 

1000 FORMAT(20G) 

«OQ.Q_.FP_R‘^Ti20F41,0) ! 

CALL FFT(N0P,NU) 

IITTJ 


5 P(L)=(TR{L)*TR(L)+TI (L)*TI (L))/(FS*F8) 

I FLCIy E 2 « i ».Q^ ^ 

wRITE(6,5000-) (PCL)ytsl,l0l2) 

LOO CL F gRMLUML. L) ^ 

call SERCH(M»^,ND2,I,PT) . 

•vRITFdOl ,5000) NM 


OUTPUT PT 


2 j=i, lui 

iicfc^iucij csUi JJ/ Jsii.'.Yj 

6 CCMI'.'.c 
L-=FS/ .r 
>^ = 2.*3. 

C n0.o£-33) •‘•F/2, 

C " V / 

CuT ? xjL..;)£^I>v,ljY 

N * 2 = N / 2 
. . f,Y2 = r.Y/2 

DO 10 I = l,*,x 

10 x(IJs(I-:JX2-n*DX 

DO 11 Isl,NY 

U...- I.U. J.a-a-J^Y,2-..U *J1Y 

OUTPUT ' ' 

PAUSE-.- hIT RETUR'^i TU GO' 

call IMTT(9b00) 

CALL Cnl.NDOCO. , 760,#0.,760. J 

call TrtirMDC(0»7‘iO,0#7flO) 

CALL-C.QlLmj?-C-X> Y,£^JjX.»-5 

CALL FIMTT(0,780) 

20 0JQ_F OR^Ji I C S DPDI N T_ I s 5 « 1 6 ) . 

3000 F0HMAT(SSTR£AK MIDPOINT ISS,I6) 

^-COJiUNUE 

END 

-S LLsf.HiTiN£ FF_Ti\i^;, ULL 

CO^'^CN/FTR/ XR (202^ XI (202^1) 

N2 = N/2 .. 

NJl =NU-1 

DO 100 L=T,NU 

102 DO IQ l 1=1. N2 

PsISITR(K/2<^*NUl ,NU) 

APG = o.2e3ie5*P/FL0.AT(NJ 

C=C0SCARG) 

S = SIMARG.J 

KlsK+1 


TRsXR(KlN2)*C+XI(KlN2)*S 

1 1 = JO. ( KJ N2 ) *C- XR { til_N21*S 

XR(K1N2)=XR (Kl )-TR 

X IC*aty2 ) =XI (Kl) -LI 

XR(K1)SXR(K1)+TR 


m 


•VOlEKillliSnli 





1 = 1 = 1 TS Cr.-2 , +1 

IMI.LE.k) go to 103 



TI = i(I (K) 

K - { O = ^ c ( I ) 

* I ( ►• ) = X 1 ( I) 

_ - Xt(I)=TR 

* I { I) = T I 

ini noniiriuL 

n’t T JRN 
E’-‘D 

Pur.CTION IBITe(J,\U) 
ieiTR=0 

DQ-2a.Q-.,I «.l Wil. ) 

J2=Jl/2 

I>=lTR = lBIT.R*2t ( J1-2*J2) 

200 J1=J2 

. _ . RET Or.'j . . ... 

£•■.0 

S 'Ja?-D l! T i:^^£P.C tiUL»rl4-\D2 » J ^ JEU 

CnwC-g/SEP/PCl 012) 

PT=,5£-&9 

DO 1 1=1,1012 

Pl=PI.N02-n 

IF(Pl.GT.PT) PTsPl fNws.\D2-I 

1 — 

FTTs.SE-OP 

. I'^(PT.LE.ptT) .%PITE(101,P) j. .. .. 

9 FOP^ATUnO max F0U%Di,2l6) 

RE TUP', 

E^O 

su rroutinf cq htojr ( x.y. z,nxx.nyy,n; 

c 

C IMPLICIT DOUBLE PRECISION (A-H^O-Z) 

REAL 

. CONTOU^ PLOrCF Z=>Tx,V) 


DIMENSION Ci(5Uj.INi5UL/LAB_EL(5LIi 

* XP(4j),YP(«),XV’(5)f YV(5),ZV(5) 

CQF'MQN /PLTCQM/ NPLTF, IPLTC,1PFLG .1PPEN _ 




CCONTAB(30),INTNSTY(30)»LABL(30),ICONOPS(5) 

_PAJLA_JL ABL OJ^ 1 * 1^20 2 H 3 _#. 2 H« 2 H 5 , 2 h^ 

* 2H9 ,2h10,2H11,2h12,2HJ3,2H1«,2H15,2H16,2H17,2h18,2H19,2H20/ 


OPTIONS IN ICONOPS(I) 

'4PUTE CONTOUR TABL 


1 * ABSOLUTE VALUE SPECIFIED IN CONTAB 


p o o o o o ri 

I 


= -£lativL valjE s‘e:kiel) r c:-ta= 

= COMPUTE INTfc\SITIE5 

= . I\TE-«SIT I.. V-Lw£S 5-ECIFlLJ I. J\T',£T> 
= 03 \0T PkJi.T CC’^TU’..-^ TAeUE» 

= PPI',T COMTUR TA 5 LE 

= scales A'.O pGpl,EP 

= INCLUDE BORDER 
= I .CLU3E scale DIvI 3I'.-\S 
= - I’.CLUDE-^CALE f.-J^'oErS 
= C-flT CCMOUR LABELS 
s CQMOL'R LABELS SPECIFIED 
= SETUP COfJTOUR LABELS 


DATA ICSIZ2, IC5IZ3 /3H 01, 3h 02/ 

..DATil XRMX, Y£.V.X-/5D , , 5 0 . , 75 0 . , 7S0 , 

DATA IXB0,IYB0,IWIN, I^‘AX,NCVS, ICO'jGPS 
*. .73,3, 12. , Ifc. ,20, 1,0, 0,3,0/ 


C'JTPUT fXX,'JYY 

Qi_L£i_J^.,2i2 

C0MAB(n=5.*I 

= 

ERROR CHECKS 

_ \X = .NXX 

r,Y = -JYY 

IE ( NX .LT. 2 ) GO TO R5 0 

IF C f.'Y .LT. 2 ) go to P50 

. IF ( NZDI^' .LT. NX J GO TO R50 
IF ( XP‘-*X ,GT, 1 700, ) XP»‘X = 1700, 

IE ( YPwx ,GT. 1000. .1 YPHXslOOO. 

U2 = \Z^S + 1 

IF ( \ I ,GT . 50 1 N ZsSO 

IE ( NZ .LT. 2 ) NZ= 2 

_ NZ2 =L NZ/2 + I .. 

NZA = \Z/a 

' COMPUT£’~MIN'*AN;b max’ 

XMTN = xm 

XMAX s X(NX) 

ymjn = Y(i; 

YMAX = Y{NY) 

_ ZMjN ^L.E3Q . 

Z^AX = 

' DO~T30 1 = 1. NX 

DO 130 Jsl,NY 

■ Z_MJN S_M INI L_ZMlli,_2.( 1 j JLJ 

ZMAX s OMAXl ( ZMAX, Z(I,J) ) 

DL AB c ( f f X P MX» XPMM)/(NX-1 ) ) **2 

« ♦ ( (YPMx-YPMS)/(Ny"-1 ) )**2 ) * .5 


contour option SETUP 



r 


r 

140 

I = ico\:°scn 

IP ( I ,NE. 1 ) GO TO 15a 
1-.4 i = d;.z - ... 

■■ 

1 U4 

CV(I) = CuM4d(I) 

C.r.Z) = 2.*C.(\Z-1} - CV(r.Z-2) 

I = 1 

I = 1 +. 1 . 

!► ( Cv(I-l) .GT. CV(l) ) GO TO 955 
IE £_ 1 *. LI NC.y3 I GO TO 145 


145 

1 

150 

GO TO 170 

IP ( I .GE. 2 ) GO TO 160 
DZ = (Z'^AX-Z^IM/l 00. 

. DO. 155 J = d^2-. 

V— 

155 

CV'(I) s CO\TAB(I)*OZ + Z'MN 

T s 1 


160 

GO TO 145 

_ C ..s fJZ - 1 ' 

DZ = (Z^AX-Zvifo/C 

1 ' = INT tLQGlOfDZl+l OCO- 1-1001 . 

> 


C = 10.**I 

P7 r T\T C DZ/C ) * C 

) 


ZP^N = INT ( Z*'‘IN/DZ ) • DZ 

IF LZI^IN ^LT. .0. .1 ZPN«\ = ZPMN-DZ 

r~ C = ZPVrg 

Cvm r ZPVN 

' DO 165 I = W-NZ 

C g C t DZ _ 

_ 

165 

170... 

cy(i+i)= c 

ZF. =. 100../(CV(NZ)-CV(U) . 

Z5 = - CV(1)*ZP 

DO 172 I=1/NZ 

i; 

172 

LABEL(I) = LA3L(I) 

I r ICONOPSfB) 

1 , 

~ i tT' 

IF ( I ,NE, 2 ) GO TO 160 
I ._ a.O . . . _ 

I =1 + 1 

C = CV(I)*ZP + ZB . 

ENCODE ( u, 9bH, LABELd) ) C 

^ IP ( I .LT. NZ ) GO TO 174 


160 

I = IC0NOPS(2) 

IP ( I .NE. 0 1 GO TO 190 

ZP = (Imax-Imin)/(CV(nZ)-CV(D) 

ZB K IMIN - ZF*CV(1) 



DO 164 1=1 »NZ 


1 84 

iNdD = "ZB + zF*cvrn ■ ■ ■ - - -• 

f 

1 86 

DO 166 1=1, NZ, 5 

INCH s IN(I) +3 . . 

i 

190 

GO TO 195 

DO 192 1*1. NZ 

\ 


K • = INTNSTr(I) 

IF ( k.GT.28 .or. K.LT.2 ) K*16 


1 92 

I N Cl ) = K 

. — 




143 



195 i = ICC, .-5(33 

IP ( I ,E'0. 0 ) GO TO 200 

^ F!iI = tCAi:£LaJ.C.m,I = l/',Zj 

C 

' c CC\TCC“ G£‘.E--TICO 

_ 2 00 i-p = ( 

' »- =XP‘\-Yr*XSIf, 

I c = .. Y P \ - _or E * Y .»-_l fi. 

IX r 1 

- X2 s xCIXl*XF + XE 

1x6 = - I X 3 0 

^ I 5 ICO.\CPS(5) 

IP ( I .Erj, 0 ) IXB = 20 

■ 

_ 205 Y( I) = y(I)*yF + YB 

C LOOP ON X-VALUES 

210 IF ( IX. .EG. la ) -GO TO 270 

C C-i-ITPuT X2 


IX = I X <*• 1 

_ >;2 =_-XClXI*.XF t XB 

XV( 1 ) = XI 

XV (2) = XI _ 

X V' ( 3 ) = X 2 

' XVX^L) =_X2 

xv(5) s XI 

IV6 -lYBO 

1X3 = IXB ♦ 1 

... IP ( IXB .,EO.. 1 ). IXB5-IX60 
I Y si 

12 -=_i.LXTJ 

ZV(2) s Z(IX-1,IY) 

ZV(3) .= Z(IX . ,IY)__„ 

C 

. C LOUP OiM Y,-V'ALUES- 


220 IF ( lY ,Ea. NY ) GO TO 210 
YJ = Y2 


1 

I YB = 

IF .C_IYB. 

lY = 

- Y2. = 

IY8 + 1 

.GT^ 0_J. I.YBs-IYBO . . _ 

lY + 1 

YCIYJ 


YV(1) 

s 

Yl 


YVC2) 

s 

Y2 


YV(3) 

s 

Y2 


• YV(4) 

3 

Yl 


YVC5) 

5 

Yl 


ZV(4) 

S 

ZV(3) 

* _ 

ZV(5) 

3 

ZV(2) 


ZVd) 

3 

ZV(2) 

\ 

ZV(2) 

1 

Z(IX-1,IY) 


i ’ 
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z V ( 3 ) = z ( I X , m 

:• = “T 1 ( ZV(1),Z.(2),Z»(3),2.(‘J) ) 

z*. i._.(_zwn,z. czj,Zv (3},zv 3 

= FO- CC- T0L"^S £ET.>£E*, Z*-‘ i‘.n 

IC = 1 

I = .iI2 

( CV(I) ,lt. 2"' ) IC = I 

I s-ic -> :^2a 

!^ ( CV(I) ,lt, z*-' ) 1C = I 
IC = 1C - 1 - 


-23C'_ ,IC. 
C 


= -IC + i 

= CV(IC) 


IF ( c .GT. ZX ) GO TO 220 
IMEFPOLATE fop CO’vTOUR intersections 

U = 1 

3 1 £_I .t.-l 

IF ( (C-ZV( I-in*(ZVCI)-C) ,LE. 0. ) GO TO 2u0 

XP(N1 .-=_XV(I-1J 

YP(N) = YVCI-U ♦ (YV(I)-YV(I-l ) ) 

_* *_tc-2vciruj/.czym-zvii-i)) 

0 = N ♦ 1 

I I =._L_t_i 

IF ( (C-ZY(I-1 ))*CZV(I)-C) ,L£. 0. ) GO TO 2U5 

_ YP(N) = YVCI-1) . 

X®(N) = XV(I-l) 4 (XV(I)-XV(I-n ) 

. * .(c-zvci- 1 ) j/(zvm-zv(i-n ) 

h. Z N 4 I 

> IF ( I . EQ. .1 ) .GO TO 235 

IF ( N ,NE, 3 ) GO TO 230 


C PLOT THE CONTOUR SEG»^ENT 

_.25J} If_.(_IXR_^EO,. 0 J GO TO 260 _ . 

C CALL UPPEN 

__255 CALL MQVEA ( XPfn>YP(l) ) 

C 

C_'^*_NCTE: the variable INTENSIT Y F EATURE 

C ** HAS BEEN DISABLED. 

C 

C CALL D^NPEN (IN(IC)) 


C CALL DwNPEN ( IPPEN ) 

CAL L DRAWA ( X P(2), YPC 

GO TO 230 

C. 

C ANNOTATE the CONTOUR 

2b0 IF ( lYB .NE. 0 ) GO TO 
lYB s -IY60 


1A5 




C CtLL 

C-'-LL •■•UVEA ( x2-?5., Y?-30. ) 

CALL. =3Cr {vi5<ICS123) 

295 xl s^l+Dx 

IF ( XI .LT. **^Ax ) GG TO 2bD 

»2 = yp ♦ fcO, 

I = tZ-OFlFO^J 
C CALL 

C.LLL..:^L«£A £ .a2^..J2?30, _J 

CALL pBCh (I,3,ICSIZ2) 

ENCCDE (3,Q8b,I) IC 
C CALL •jPA’EN 

call mOvEA ( X2+25.4 Y2-15. ) 

CALL PbCH (I,3,ICSIZ2) 

C 

DV = ( V'^AX-Y^IfO/CYP^'X-YP*' N)*100, 

1C = I\'T (LOGIC IUY3 + 999,PQ?999'») - 1000, 

C = 10.**IC 

_ - DY = IfiT( J3r/Ct,9999P9P9 ) * C 

Y1 = IijT(YVIN/DYJ*DY + OY 

i£_L-3_-aii_j_i,_ii._i n_sjfi-D.y 

XI S xPMN - 10, 

X2 =_XP.!^ ... 

Z1 = XPVN -f XPVX - XI 

22 = XPVN ■» XPM.X « X2 . 

Y2 = YMAX 

C C 

268 C = .1 * C 

IC = IC .♦ I 

IF ( A8S(Y1*C) ,,GE, 10, ) GO TO 288 

IF (_ABS(Y2*C) .GE,_10. ) GO TO 288 . 

290 Y2 = YF*YI * VB 

C CALL U PP EZJ 

CALL y.UVEA ( XI, Y2 ) 

C call Oanpen ( IPPEN ). 

CALL DPArtA ( X2» Y2 ) 

C call UPPEN .... - 

call mOvEA ( 21, Y2 ) 

C CALL DrtNPFN f IPPFK 1 

CALL DRAWA ( 22, Y2 ) 

lF-_Lfl_.fJE.. 3 J GQ ,T0..295_ . 

S = Yi*C 

ENC0DE_L9,_9M, s 

C CALL UPPEN 

■ - --- CALL-^OVEA — r-X2-70 , -.-— Y2-5,--rr-.:: -.. 

call P8CH (V,5,ICSIZ3) 

2.9.5_ Y1 s Y 1 _t. D Y 

IF ( Y1 ,LT, YMAX ) GO TO 290 

Y2 = ..Y.2 5 0 . 

I * ezaoFiFoao 

C call UPPEN 

call MOVEA ( X2-55,, Y2 ) 



L 


\ 

i 

l 

I 

r 


« 


t 


CilL (1,3, 1CSIZ3) 

E^Cull (3,Pre,n 1C 

_.__C LALL _'r£l. 

C^lL “‘0;E4 ( 1^2-30., Y2*15. ) 

C^LL “r:*' (I,3,IC5IZ£) 

C 

.. C • 

3 C C V r r 1 . / V r 

Da..3 LiS J?1 , \Y - - 

305 v(l) = (Y(I)-Y5)*YF 
C 

RE TU-*.' 

C ‘•‘ESSiGE 

OF3 c = TM — - 

GO TO 

955 PRIM 955, ( CDMABU ) , 1 = 1 # M) 

960 STOP 

C - - .- 

963 FOR*-'iT ( 1 HI //Tl 1 0,20^SY*^50L CONTOUR VALUE 

* J / ,.(TI1 Q , kZu.EA3 •BJ. 1 

«?54J FOR '‘AT ( Fa, I ) 

965 FOPvat (__f5.2 ) . 

956 FCR--AT ( 13 ) 

957 FCP.-^AT ( IHl ,2DX22^'E9PUP JN_.CCf.TOUH USAGE / 

* 1 Ox3-‘.x=, la , 1 0X3^-».Y=, la, ICxSh'-'DIwz, !«// ) 

0 fi 5 F Q P 1 ‘ AT C 2 0 .x_ 3 2 H £ 0^ 1C. y ®_ ,V A L U E S _N Q T . *•' 0 N 0 T C I C A L L Y 

• li-. INCREASING /7 c5x,5E20.6) ) 

End 

S“=®C.TrjE P.cC'^' IL/N',*') 

C 

CI''ENsTcN”LTao) ,K’U3) , IDCRiTTrrfO 

DATA ID /IHJ , l^.:2rlH3/!Ha,lH5,IH6,lH7,lHe,lh9/ 
DATA lA /IhO, IMI , 1m2, 1H3/ 

C 

CALL AivVODE 

L£_X L £ -J2 1_£ LI U£Ii 

IF (N.GT.iJ) NNSU 

GO TC (10,20, 30, ... 

10 WRITE (102,15) L(l) 

_ 15 format (An 

RETURN 

~ 20 i-.PTTF -M0?.?5) Lfn 

25 FOR'^AT (A2) 

RETURN.:.. .. . 

30 WRITE (102,35) L(l) 

_ J5 fUJi.lAl (A3) 

return 

gp Y.P1TE (n Z iqSJ.iUJ 

a5 format (Aa) 


return 

END 


I 
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11. Task 10: Single Aerosol Response 


The design and fabrication of an apparatus for generating 
and propelling single aerosols through the lidar beam in 
order to determine the lidar response has been previously 
described. The technique for measuring this response 
requires only that the lidar response to a particle crossing 
a region of the beam of known relative calibration be 
determined. When this is known, each relatively determined 
contour in the focal volume map can be absolutely 
determined. 

This response can be ascertained from an ensemble of time 
domain single particle responses obtained from particles 
near the beam response maximum. A display mode which simply 
accumulates these responses in superposition can then be 
read to determine the maximum response, from a group of 
identically backscattering particles. Care must be taken to 
assure that only single particle returns are used. If in 
doubt, a histogram display of pulse heights should indicate 
any multiple responses. 

The state of calibration of the lidar system must be known 
during these measurements in order that, under different 
calibrations, the single particle response can be determined 
from the measured value by scaling. That is, when the 
system is recalibrated with a waist measurement, the new 
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single particle response is just the ratio of the waist 
calibrations times the old single particle response. 

Single particle signals from the single particle generator 
have been obtained and are seen in Figure 10.1 showing a 
spectrum analyzer trace. Heights of these pearcs are 
proportional to the particle scattering cross-section and 
are also dependent on where the particle crosses the focal 
volume. If all particles have the same backscatter cross- 
section, the highest peaks indicate the single particle 
response of the system as each particle crosses the center 
of the beam waist. 

Single Particle Algorithm: An algorithm for measurement of 

Beta under conditions such that only single particles give 
significant scattering from the focal volume has been 
developed. It can be shown that, utilizing ordinary volume 
mode processing techniques, much smaller backscatter values 
can be determined in the single particle mode. Under the 
assumption that the response of the system at each point of 
the focal volume is known, the following equation yields 
Beta per unit flight path from the distribution of peak 
signals collected over some time interval (and therefore 
over some spatial region): 
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Here 


N = number of particles seen 
0^ = backscatter cross-section 

= probability of o^. 

a^. = focal area sensitive to o. 

The are found by calibration and mapping of the focal 

volume, while are determined by fitting the following 

equation for the peak signal distribution, Sj, to the 
measured peak signal distribution using a least squares fit 


A.. = Cross-sectional area of focal volume 
sensitive to and yielding a peak 
signal S. 

J 

n . = N <u./Z u).a . 

1 1 V J 

c number of particles with cross-section 
per unit area. 


The distribution 


is assumed to be parameterized by 


0 )- 

1 


B/o^ 


b 


Hence the single particle algorithm 
152 


where 0 < B < 3. 



1) assumes N determined by real time measurement 

2) assumes a ^ and A are determined by system 
calibration 

3) loops over the parameter B until Sj fits the measured 
values in a least square sense 

il) takes the successful B parameter and calculates p. 

Simulated data from a system with an ideal Gaussian beam 
focussed at 5 meters with a 10 centimeter diameter mirror 
were used to test the algorithm. Random parameters in 
the simulation were the actual distribution, and the 
particle positions. 

Table 10.1 shows the effect of the number of integration 
steps on the algorithm for three different simulated power 
law densities (B values). It is generally noted that as few 
as 5 integration steps gives reasonable results. 

Table 10.2 shows the effect of "signals seen" on the Beta 
algorithm. These results show that no noticeable 
deterioration is evident until the # of particles seen drops 
to between 10,000 and 5,000. 

Tables 10.3 and 10.4 show the effect of varying the # of 
signal and sigma bins. The general conclusion reached here 
is that the important variable is the "# of signals over the 
threshold." When this number becomes a large fraction of 
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the total of signals then the results may not be 

realistic. 

Procedure Required to Run the Single Particle Inversion 

Algorithm: 

In order to run the ys inversion algorithm the following 

steps should be followed: 

1) establish correspondence between signal bin number and 
linear signal value 

2) establish noise value, since algorithm uses s/n 

3) determine s/n threshold and set SIGL = threshold s/n 

4) set number of signal bins (M) and number of sigma bins 
(MS) - in general, number of sigma bins should be larger 
than number of signal bins - say M = 4 and MS = 7 

5) set SIGH = highest linear S/N value 

6) input a signal histogram that corresponds to the number 
of bins desired and the size of each bin in array NP; 
also as an N + 1 element include the number of particles 
above the maximum s/n. 

Other parameters of interest are: 

KI - number of integration steps 
ETA - system efficiency 
BW - bandwidth 154 



FSI - sigma bin size relative to preceding bin size, 
i.e., FSI = 2 means that each bin doubles in 
size from the first bin interval to the last 
bin interval. 

The algorithm as currently coded uses a theoretical Area 
versus s/o table. The experimental determined Area versus 
a/o' curve can be used to replace this "Area table.” 
Subroutine AREA does the theoretical calculation and should 
be removed. Array AR (101) is used to store the "Area 
table." 
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Table 10.1 


Effect of 

r- Steps 
30 
20 
10 
5 

30 

20 

10 

5 

30 

20 

10 

5 


- of Integration steps 

h (XIO'^) 
2.459949 
2.459994 
2.460194 
2.461007 

2.180584 
2.180605 
2.180709 
2.181125 

2.108201 
2.108227 
2.108355 
2.108887 


Beta algoritiii 

B Value 

2.5 

2.5 

2.5 

2.5 

1.5 

1.5 

1.5 

1.5 

.5 

.5 

.5 

.5 


Fixed Parameters 

Particle Site (p) - (1* " 1.316228) 
Cross section (geometric) -(n-10ii) 
Focal length - 5.0 

# Sigma bins - 7 
f Signal bins - 3 

# Signals seen - 25,000 
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Table 10.2 


Effect off signals seen" on (i algorithm 


# signals Selected B (tr’./^signal s) X 10 


30,000 

1.5 

B. 7223366 

25,000 

1.5 

8.722336 

20,000 

1.5 • 

8.722335 

15,000 

1.5 

8.72234 

10,000 

1.5 

8.722335 

5,000 

1.4 

8.663662 

2,000 

1.3 

8.613045 

1,000 

1.3 

8.613041 

500 

1.5 

8.72233 


Fixed parameters 

No. of integration steps 
Particle size (y) 

Cross section (geometric) 
Simulated B 
Focal length 
#sigma bins 
#signal bins 


- 30 

- (1 - 3.16228) 
“(rr -10 t) 

- 1.5 

- 5.0 

- 7 

- 3 
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Table 10.3 


Effect of varying ^ signal bins, MS = 10 (# sigma bins) 


t sianal bins 

E (X 10'^) 

Selected B 

^ sianal s over threshold 

2 

2.197394 

1.6 

24,218 

3 

2.197394 

1.6 

23,237 

4 

2.180583 

1.5 

21,538 

5 

2.180583 

1.5 

18,554 

6 

2.180583 

1.5 

14,325 

7 

2.180583 

1.5 

9,057 

8 

2.180583 

1.5 

4,224 

9 

2.180583 

1.5 

1,011 

10 

2.180583 

1.5 

1,011 

Total r Signals 

- 25,000 



B 

- 1.5 



# Integration steps 

- 30 



Particle size (p) 

- (1 - 3.16228) 



Cross section (geometric) -(it -10 -n) 



Focal length 

- 5.0 




# 
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Table 10.4 


Effect of varying #' of sigma bins 
M=3 signal bins) 


a sigma bins 

6 (X 10"^) 

Selected B 

# signals over threshol 

10 

2.197394 

1.6 

23287 

9 

2.180582 

1.5 

21744 

8 

2.180585 

1.5 

19054 

7 

2.180584 

1.5 

15058 

6 

2.180592 

1.5 

9839 

5 

2.180606 

1.5 

4830 

4 

2.180649 

1.5 

1243 

Total # signals 

- 25,000 



B 

- 1.5 



a integration steps 

- 30 



Particle size (y) 

- (1.-3.16228) 



Cross section (geometric) - (ti -10 r) 



Focal length 

- 5.0 
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12. Task 11: Software for Transfer to Sigma V Computer 

Software was developed to transfer backscatter data which 
was' collected and stored on disk files on the Beta signal 
processor. The LSI-11 processor was interfaced with the 
Sigma V Computer and the following functions automatically 
performed: sign-on, opening of the data file on the Sigma, 

transfer of data, close data file, and sign-off. The 
transfer of more than one data disc is also possible. 
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13. Task 12: Interface of HP 1 POOL Computer with 


Peripheral Devices 

API has developed necessary software to interface a Hewlett- 
Packard 1000L computer with peripheral devices for the 
purpose of transferring, processing, displaying, and storage 
of recorded data. Three routines have been written for the 
project . 

• A listing routine to provide Source code listings 
using the system plotter. This routine is fully 
operational. A list of this code follows. 

• A routine to real LeCroy Digitizer data and transfer 
to the hard disk. 

• A routine to read data from the hard disk and make 
plots. A digitizer routine is designed to 
facilitate data collection in the field. Data is 
transferred from the digitizer in a block whose 
starting point and length is within the total sample 
frame and is user defined for each frame. This may 
be as much or as little as desired. The data is 
packed two data points per computer memory word and 
transferred to hard disk on a track and sector 
basis. This design bypasses the operating systems 
file management but is an order of magnitude faster 
with a typical data transfer taking 4 to 6 seconds. 
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The disk read and plot routine is intended only for 
data verification and as an example of how the user 
should read data from the disk to any processing 
program which becomes necessary. 
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1^. Conclusion 


For this contract Applied Research, Inc. has accomplished a 
broad range of tasks in support of the NASA/MSFC Laser 
Doppler Velocimeter project. These tasks involved hardware 
procurement and fabrication for the laser system, software 
development, system calibration, flight support, and data 
analysis. A two color LDV system design analysis was 
accomplished and an adaptative filtering signal processing 
technique was evaluated, both in support of future systems. 
Notable in the area of hardware development were the design 
and fabrication of a single particle generator and a 
spinning wire target, both used for coherent focal volume 
calibration. Coherent focal volume mapping was accomplished 
for the first time on LDV systems. A mathematically 
formulated defination of the relationship between hard 
target calibration and atmospheric return was formulated in 
Lg£.j|.. Many of these efforts serve, to enhance the validity 
of backscatter measurements in the volume mode (high return 
case). New in this effort was the implementation of the 
single particle mode of backscatter measurement for the case 
of low backscatter, which, in addition, yields the 
backscatter cross-section distribution. This mode may be 
capable of providing backscatter measurements when the 
volume mode is not sufficiently sensitive. 







